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Abstract 



We consider spatial coarse-graining in statistical ensembles of non-selfintersecting and one- 
fold selfintersecting center- vortex loops as they emerge in the confining phase of SU(2) 
Yang-Mills thermodynamics. This coarse-graining is due to a noisy environment and 
described by a curve shrinking flow of center-vortex loops locally embedded in a two- 
dimensional flat plane. The renormalization-group flow of an effective 'action', which is 
defined in purely geometric terms, is driven by the curve shrinking evolution. 

In the case of non-selfintersecting center-vortex loops, we observe critical behavior of 
the effective 'action' as soon as the center-vortex loops vanish from the spectrum of the 
confining phase due to curve shrinking. This suggest the existence of an asymptotic mass 
gap. 

An entirely unexpected behavior in the ensemble of one-fold selfintersecting center- 
vortex loops is connected with the spontaneous emergence of order. We speculate that 
the physics of planar, one-fold selfintersecting center-vortex loops to be relevant for two- 
dimensional systems exhibiting high-temperature superconductivity. 



Zusammenfassung 

Die Anregungen der konfinierten Phase in der thermodynamischen Behandlung der SU(2) 
Yang-Mills Theorie sind Zentrumsvortexschlaufen welche aufgrund der Wechselwirkung mit 
einer rauschenden Umgebung Schrumpfungsprozess unterliegen. Wir betrachten statisti- 
sche Ensemble von Zentrumsvortexschlaufen ohne und mit einfachen Schnittpunkt welche 
in einer flachen zweidimensionalen Ebene lokal eingebettet sind. Der Schrumpfungspro- 
zess von eingebetteten Zentrumsvortexschlaufen wird durch eine Diffusionsgleichung be- 
schrieben. Der Renormierungsgruppenfluss einer in rein geometrischen Grofien definierten 
effektiven „Wirkung" wird durch die Evolution schrumpfender Kurven bestimmt. 

Im Falle von Zentrumsvortexschlaufen ohne Schnittpunkt beobachten wir ein kritisches 
Verhalten der effektiven „Wirkung" sowie die Vortexschlaufen aufgrund des Schrumpfungs- 
prozesses aus dem Spektrum der konfinierten Phase verschwinden. Dies legt die Existenz 
eines asymptotischen Massen-Gaps nahe. 

Ein vollkommen unerwartetes Verhalten im Ensemble von Zentrumsvortexschlaufen 
mit einfachem Schnittpunkt steht in engem Zusammenhang mit dem spontanen Auftreten 
von Ordnung. Wir vermuten, dass die Physik ebener Zentrumsvortexschlaufen relevant ist 
fur die Beschreibung zweidimensionaler Systeme, welche die Eigenschaft der Hochtempe- 
ratursupraleitung aufweisen. 
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Chapter 1 

Introduction 



The importance of Yang-Mills theories in mathematical and theoretical physics is generally 
acknowledged. Yang-Mills gauge theories are the cornerstone of quantum field theories in 
the Standard Model of Particle Physics: Besides gravity, all fundamental interactions are 
incorporated as gauge symmetries in the Standard Model. Although it has been examined 
in the framework of perturbation theory due to the enormous complexity implied in the 
full story of (especially non-Abelian) gauge theories, the Standard Model has produced a 
lot of striking results and predictions. There are many examples, such as the explanation 
of the anomalous magnetic moment of the electron, the feature of asymptotic freedom of 
Quantum Chromodynamics in the high energy limit, or the prediction of flavor-changing 
neutral currents in electroweak processes [T]. However, there are still a number of unsolved 
mathematical problems and unexplained experimental observations. Among those are: 
The necessity of an asymptotic mass gap and a rigorous proof of color confinement in 
pure Yang-Mills theory [2]. In the Standard Model, the assumption of a zero rest mass 
of the neutrino is refuted by the observation of neutrino oscillations |3j and the double 
(5 decay [I]. These observations indicate a small, finite rest mass that also cannot be 
excluded by recent experiments measuring the spectrum of the single decay of tritium 
nuclei near the endpoint [51 IB]- Furthermore, the Standard Model does not provide for an 
explanation of Dark Matter and Dark Energy that account for about 96% of the energy 
density in the present universe, and the predicted Higgs particle has evaded experimental 
detection so far. Moreover, the perturbation series of four- dimensional quantum field 
theories is most likely an asymptotic series; the fact that a perturbative calculation of the 
thermodynamical pressure cannot be driven beyond order g 5 in the coupling constant due 
to the weak screening of the magnetic sector causing infrared instabilities [7], could be 
shown for Quantum Chromodynamics at finite temperature. 

Since perturbation theory is an expansion in powers of a necessarily small coupling 
constant about a trivial a priori estimate for the vacuum of the theory, it fails to describe 
strongly coupled physics as well as the according nontrivial vacuum state. This vacuum 
is certainly composed of finite-action solitonic solutions of the classical Yang-Mills action. 
The so called instantons are topologically nontrivial objects in pure Yang-Mills theory 
that describe tunneling processes between topological distinct vacua, e.g. [8]. Their weight 
possesses an essential zero at vanishing coupling, and thus instanton contributions to the 
partition function of the theory are completely ignored by perturbation theory. Instantons 
at finite temperature are called calorons. 

Therefore, we are advised to consider a nonperturbative approach to gauge theories. 
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CHAPTER 1. INTRODUCTION 



Such a treatment has already been proven successful in terms of an effective theory for 
superconductivity [9]. An analytical and nonperturbative approach to SU(2) Yang-Mills 
thermodynamics was developed in [10]. In this approach the basic idea is to subject the 
highly complex dynamics of the topologically nontrivial field configurations to a spatial 
coarse graining that leads to the emergence of macroscopic scalar fields, and pure gauges. 
Due to nontrivial (thermal) ground states, the fundamental gauge symmetry is broken 
successively as temperature decreases. As a consequence, Yang-Mills thermodynamics 
occurs in three phases: the deconfining, the preconfining and the confining phase. The 
latter, in which we are primarily interested in this thesis, exhibits three unexpected results. 
These are the exact vanishing of the energy density and the pressure of the ground state 
at zero temperature, the Hagedorn character of the preconfining-confining phase transition 
and the spin-1/2 nature of the massless and massive excitations in the confining phase. 

The ground state of the deconfining phase is composed of interacting calorons and an- 
ticalorons and exhibits negative pressure. The propagating excitations within that phase 
are two massive gauge modes - due to the dynamically broken SU(2) - and one massless 
gauge mode. As temperature decreases, the likeliness for calorons and anticalorons to dis- 
sociate into (BPS saturated) magnetic monopoles and antimonopoles increases strongly in 
the vicinity of a critical temperature. The ground state of the preconfining (or magnetic) 
phase starts to form by the pairwise condensation of monopoles and antimonopoles. Exci- 
tations in that phase are propagating dual gauge modes of mass m D (dynamically broken 
U(1)d)- Unstable defects of the magnetic ground state are closed magnetic flux lines of 
finite core size d that collapse as soon as they are created. This is because, as long as 
d > 0, the pressure inside the vortex loop is more negative than outside, thus leading to 
the contraction of the vortex loop. The magnetic phase exhibits negative pressure. At 
the Hagedorn transition towards the confining (also called center) phase, a complete de- 
coupling of the gauge fields takes place. To put it more precisely, by the decay of the 
magnetic ground state into selfintersecting and non-selfintersecting center- vortex loops the 
mass of the dual gauge field diverges and the core size of center-vortex loops vanishes, see 
also [HUE]. As a result of d — > 0, the negative pressure P is confined to the vanishing 
vortex core. This implies that center-vortex loops become stable particle-like excitations 
with P = 0. These solitons are classified according to their center charge and the number 
of selfintersections N, see Fig. 11.11 The mass of an A r -fold selfintersecting soliton is NAq, 
where Ac is the Yang-Mills scale. Topologically, solitons with non-vanishing N are stable 
in the absence of external gauge modes coupling to the charges at the intersection points. 
On the other hand, for N = 0, there is no topological reason for stability. 

Now consider a situation where a planar center-vortex loop, which is a (local) em- 
bedding of a center-vortex loop into a two-dimensional flat and spatial plane, still has 
non- vanishing core size d > and the mass of the dual gauge field m D is still finite due 
to a noisy environment which locally resolves the otherwise infinitely thin vortex. In this 
case, the pressure P is locally nonzero and the soliton starts shrinking. Such a situation 
is described by a curve shortening flow in the (dimensionless) parameter r. Here, r is a 
variable measuring the decrease of externally provided resolving power applied to the sys- 
tem. There is a functional dependence of r on the corresponding resolution Q (momentum 
transfer). For an isolated SU(2) theory the role of the environment is played by the sectors 
with N > 0. If the confining SU(2) is part of a world with additional gauge symmetries, 
then a portion of such an environment arises from the mixing between the corresponding 
gauge groups. Either way, a center-vortex loop acquires a finite core size and as a conse- 
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N = 3 



Figure 1.1: The topologies of solitonic excitations with up to N = 3 selfintersections for 
an SU(2) Yang-Mills theory in the confining phase. A magnetic monopole of charge +1 
or —1 is located at each point where center- flux lines intersect. Solitons with N = 
are unstable in presence of a noisy environment, whereas solitons with N = 1 are always 
stable. Excitations with N > 1 are unstable if subjected to mixing with theories possessing 
propagating gauge fields. 

quence, a finite mass for the N = soliton by frequent interaction with the environment 
after it was generated by a process that was subject to an inherent, finite resolution Qq. 

Knot-like structures are relevant in a number of chemical, biological and physical sys- 
tems |13| . e.g. in polymer physics, particularly in molecular biology, in type-II supercon- 
ductor, where string-like vortices confine magnetic fields to the cores of the vortex- like 
structures, in superfluid helium ( 4 He), as well as in liquid crystals. As early as 1897 Lord 
Kelvin proposed that elementary particles - at that time atoms were considered to be el- 
ementary by Kelvin and others - should be described as knotted lines of vortex tubes in 
a medium (the aether) |14j . As we know now, the point particle interpretation of Quan- 
tum Mechanics appears to be a much more elegant and efficient framework to describe 
the physics of atoms and molecules. But at the same time, the notion of an electron as a 
spinning point particle, albeit an excellent description in a bulk of physical situations in 
atoms, colliders and condensed matter systems, causes theoretical and experimental incon- 
sistencies. On the one hand, there is the problem of diverging classical self-energy of the 
electron. On the other hand, the unexpected explosive behavior in recent high-temperature 
plasma experiments [T51C16] and the strong correlations of electrons in two-dimensional pla- 
nar systems [9] are indications of non-local effects possibly related to the extended spatial 
structure of the electron. Also, recent theoretical developments revive Kelvin's description 
of elementary particles as non-local knot-like entities. In [TH [T7], the argument is that 
confining strings, tied into stable knotted solitons, exist when decomposing the gauge field 
in the low-energy domain of four- dimensional SU(2) Yang-Mills theory. 

According to the approach in [10], we tend to interpret one-fold selfintersecting center- 
vortex loops as electrons and accordingly non-selfintersecting center-vortex loops as neu- 
trinos. This implies that the Yang-Mills scale Ac must be set equal to the electron mass 
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CHAPTER 1. INTRODUCTION 




Figure 1.2: Points on the center flux lines moving oppositely on a line perpendicular to the 
bisecting line of the angle a with velocity modulus v\ . For sufficiently small a the velocity 
modulus t>2 of the intersection point is superluminal: vi = v\ cot ^. 

m e = 511 keV. The spin-1/2 nature of a center- vortex loop is a consequence of its two-fold 
degeneracy with respect to the direction of flux which is lifted in the presence of an electric 
or magnetic background field. It should be noticed for a selfintersecting center-vortex loop 
that, as long as both wings of center flux are of finite size, a spatial shift of the intersection 
point requires a negligible amount of energy only. In particular, if the inner angle a be- 
tween in- and out-going center-flux at the intersection is sufficiently small, then a motion of 
points on the vortex line that is directed perpendicularly to the bisecting line of the angle 
q easily generates a velocity of the intersection point which exceeds the speed of light, 
e.g. Fig. 11.21 Here, it should be considered that the path-integral formulation of Quantum 
Mechanics admits such superluminal motion in the sense that the according trajectories 
contribute to transition amplitudes [T8] , 

The purpose of this thesis is to treat the behavior of N = and N = 1 center- vortex 
loops under curve shrinking as a Wilsonian renormalization-group flow governed by an 
effective 'action'. The term 'action' is slightly misleading since we do not aim at a time 
evolution of the system by demanding stationarity of the 'action' under curve variation. 
We consider resolution dependent statistical ensembles in the presence of an environment 
represented by a parameter a. The corresponding weight-functional for the members of the 
ensemble, written as the exponential of an 'action', is defined in purely geometric terms. 
In turn the resolution dependence of the 'action' is determined by the curve shrinking 
flow. The 'action' possesses a natural decomposition into a conformal and a non-conformal 
factor. We consider the partition function of a given ensemble of planar curves to be 
invariant under the condition of changing the resolution. Once the evolution of the weight- 
functional is determined, we are able to compute the resolution dependence of 'observables' 
as ensemble averages of (local or non-local) operators. 

What we observe is that the N = sector becomes unresolvable from a finite resolution 
Q* downward. That is, as a consequence of a noisy environment planar N = center- vortex 
loops shrink to points with circular limiting shape within a finite decrease of resolving power 
and thus disappear from the spectrum of the confining phase of SU(2) Yang-Mills theory 
for a resolution smaller than the critical Q*. Since center- vortex loops with N > have 
finite mass this generates an asymptotic mass gap. We show that the observed transition 
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to the conformal limit of vanishing curve length is a critical phenomenon with a mean-field 
exponent of the coefficient associated with the non-conformal factor. For the N = 1 sector 
we observe the unexpected behavior that, starting from a finite value, the entropy of the 
system decreases to an almost zero value as the resolving power is lowered: the ensemble 
evolves into a highly ordered state in a sense that only a single curve survives the process 
of coarse-graining. 

The thesis is organized as follows: Chapter Ogives a brief outline of the effective theory 
of thermalized SU(2) Yang-Mills dynamics in all of its phases as it is developed in [10J. 
Chapter [3] provides prerequisites for the mathematics of curve shrinking flows in two and 
three space dimensions. Chapters [U and [5] investigate the N = sector and respectively, 
the N = 1 sector. In Sections 14.11 and 15.11 we explain our philosophy underlying the 
statistics of geometric fluctuations and how the renormalization-group flow of the effective 
'action' is driven by the curve- shrinking evolution of the members of a given ensemble of 
N = and respectively, N = 1 center-vortex loops. In Sec. 14.21 we explain our numerical 
analysis concerning the computation of the effective action and the variance of the 'center 
of mass' which is compared to Heisenberg's uncertainty relation. In Sec. 15.21 we elucidate 
our numerical analysis concerning the computation of the effective 'action', the variance of 
the location of the selfmtersection, and the resolution-dependent entropy associated with 
a given ensemble. Chapter [6] deals with electrons which are interpreted as center-vortex 
loops with one selfmtersection. In Sec. 16.11 we give reasons for this interpretation, and 
in Sec. 16.21 we consider strongly correlated systems of electrons in cuprates exhibiting 
high-temperature superconductivity and the new class of recently discovered iron-based 
high-temperature superconductors. Chapter [7] gives a short summary of our findings. 



Chapter 2 

Brief review of SU(2) Yang-Mills 
thermodynamics 

In this section, we give a short outline of the analytical and nonperturbative approach 
to SU(2) Yang-Mills thermodynamics as it is developed in [10]. The basic idea is to 
subject the highly complex dynamics of topologically nontrivial field configurations to a 
spatial coarse-graining that is described by emergent macroscopic scalar fields, one for each 
phase. Conceptually, this approach is similar to the macroscopic Landau-Ginzburg theory 
of superconductivity. Although we are only concerned with the confining phase in this 
work, we start our outline of |10| in the deconfining phase at high temperatures which 
leads us by consecutive phase transitions to the confining phase. 



2.1 Basics of thermal Yang-Mills theory 

Yang-Mills theories are non-Abelian gauge theories whose Lagrangian is demanded to be 
invariant under local gauge transformations. In this thesis, we restrict ourselves to the case 
of SU(2) gauge transformations. In pure Yang-Mills theory, only gauge field terms appear 
in the fundamental Lagrangian while matter fields are absent. Wick-rotating to Euclidean 
signature by t — ► —it and moving to finite temperature T, which corresponds to imaginary 
time compactified on a circle with circumstance (3 = ^, the gauge-invariant action is given 
by 



S = ^2 tr I dr / d 3 xF^F^, (2.1) 



1 

2? 



where g denotes the dimensionless coupling constant and tr the trace operation. It holds 
that (xi, X2, xs, X4) 6 IR 4 . The Yang-Mills field strength tensor is defined a£| 

F lu , = d ll A v -d v A ll -i[A ll ,A v ], (2.2) 

with the Lie-algebra valued gauge fields in the adjoint representation 

A, = A^, (2.3) 



lr The gauge coupling g is absorbed in the definition of the gauge fields. 
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where the generators a a are given by the Pauli matrices. The action density tr F^ V F^ V 
is invariant under local SU(2) gauge transformations 

Apix) ^> n(x)A^(x)&(x) + in(x)d^n(x), (2.4) 

where O is an element of SU(2). 

Instantons are localized finite-action classical solutions in Euclidean field theory The 
BPST (Belavin-Polyakov-Schwartz-Tyupkin) instanton is an (anti)selfdual, that is BPS 
(Bogomornyi-Prasad-Sommerfield) saturated, configuration solving the Euler-Lagrange 
equations D^F^ = subject to the action (|2.1I) [19] . For the covariant derivative Z? M 
of the field (ft in the adjoint representation we have 

D^ = d^-i[A^4>\. (2.5) 

The (anti)selfduality condition reads 

Ffiv = ±F liV , (2.6) 

where the dual field strength is defined as F^ = ^e^^-F^, e^ VK x being the total anti- 
symmetric tensor with £1234 = 1. (Anti)selfdual configurations saturate the BPS bound on 
the action which therefore is minimal (in a given topological sector) and of value 

S=^-\Q\, (2.7) 
where the Pontryagin index Q is a topological invariant (charge) and defined as 

Q = 3^2 J dr J d 3 x F; v F« v . (2.8) 

BPS saturated field configurations A^ have vanishing energy-momentum tensor. 

(Anti)calorons are BPS saturated, periodic-in-r configurations at finite temperature 
with finite action and topological charge Q = ±1. They are classified according to the 
eigenvalues of their Polyakov loop (time-like Wilson loop evaluated in periodic gauge) at 
spatial infinity. An (anti)caloron is said to be of trivial holonomy, if its Polyakov loop, eval- 
uated at spatial infinity, is an element of the center of the gauge group. Otherwise it is said 
to have nontrivial holonomy. The Harrington- Shepard (HS) (anti)caloron is a periodic-in-r 
instanton in singular gauge with topological charge Q = ±1 and trivial holonomy, whereas 
the Lee-Lu-Kraan-van Baal (LLKvB) (anti)caloron is of nontrivial holonomy. Descriptively, 
trivial holonomy means that the caloron has no substructure. The LLKvB (anti)caloron 
contains BPS magnetic monopoles constituents which, by virtue of quantum corrections 
|20J, are subject to an attractive interaction in the case of small holonomy and to a repul- 
sive interaction for large holonomy. In the case of large holonomy, the repulsion leads to 
a dissociation of the caloron into a pair of a screened magnetic monopole and antimono- 
pole. On the other hand, for small holonomy, the (anti)caloron collapses back to the stable 
configuration of a HS (anti)caloron by annihilation of their BPS monopole constituents. 
Thus, single LLKvB calorons are unstable under quantum deformations. 



2.2. THE DECONFINING PHASE 
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2.2 The deconfining phase 

The complex microscopic dynamics in Yang-Mills theory does not seem to allow for a 
direct analytic calculation of macroscopic quantities in terms of the fundamental gauge 
fields. A spatial coarse-graining, that is the computation of a spatial average over the 
sector of topologically nontrivial, BPS saturated field configurations of trivial holonomy 
turns out to be a feasible and thermodynamically exhaustive approach. The coarse-graining 
procedure is described in terms of a macroscopic adjoint field cp. In order to characterize 
the macroscopic ground state, cp has to satisfy for following conditions: 

(i) due to spatial isotropy and homogeneity cp must be a Lorentz scalar; 

(ii) homogeneity of the ground state implies that the modulus of cp is independent of space 
and time. A dynamically generated Yang-Mills scale A enters this modulus as a parameter; 

(iii) cp is a composite of local fields and therefore has to transform under the adjoint rep- 
resentation, because in pure Yang-Mills theory all local fields transform under the adjoint 
representation of the gauge group; 

(iv) only the color orientation of cp in a given gauge, also referred to cp's phase, depends 
on r. Since cp is constructed from (anti)calorons, which are periodic in Euclidean time, c/>'s 
phase is also periodic in r, and since the classical caloron action S = ^ is independent of 
temperature cp's phase is not explicitly time dependent. The computation of the phase of 
cp does not require any information about the Yang-Mills scale. 

Consequently, the field can be written as 

In [10j[2TJ[22] equations of motion for the phase and modulus of the spatially homogeneous, 
composite, emergent adjoint scalar field cp obeying the above conditions are derived. The 
(non-perturbatively) temperature dependent modulus is given by 



A 3 ^ / A| 



The corresponding action is found to be 

S4> = taJ dr J d 3 x ((^) 2 + A|(/>- 2 ), (2.11) 

where cp^ 1 = r^. The field cp turns out to be quantum mechanically and statistically 

inerid. It serves as a spatially homogeneous background for the topologically trivial (Q = 0) 
sector of the coarse-grained, propagating gauge fields a M . In Eq. (|2.11h interactions between 
calorons are not yet included. This is done via minimal coupling by substituting 

V- > - D M^ = V + M^OmI- ( 2 - 12 ) 

The interactions are mediated by the topologically trivial fields that change the holonomy 
of the (anti)calorons and subsequently induce interactions between the magnetic monopole 



2 This can be checked by direct computation but also is implied by the fact that a spatial average over 
non-propagating gauge fields must generate a composite that itself is not propagating. 
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constituents of nontrivial holonomy (anti)calorons. The action for the minimally coupled 
fields is given by 

S = tr j\r f d 3 x G^ + {D^f + A 6 0" 2 ), (2.13) 

where the field strength is G^ v = ^-(d^a® — d v a a ^ — ef abc a l \afy and e denotes the effective 
gauge coupling which determines the strength of interaction between topologically trivial 
gauge field fluctuations and the macroscopic field (f>. Due to Lorentz invariance, gauge 
invariance, perturbative renormalization, and the inertness of (f> the action (j2 .13(1 is unique. 
The topologically trivial sector is written as a decomposition 

a„ = aff + 6a„, (2.14) 

where aff is a pure-gauge solution of the equations of motion for following from action 
(I2.13p . and <5a M is a (periodic) finite-curvature propagating fluctuation. The pressure P E S 
and energy density p 9 ^ of the ground state, following from the energy-momentum tensor, 
read 

P°E = -Pe = -4vrA|T. (2.15) 

Microscopically, the negative ground state pressure arises from the creation and annihi- 
lation of BPS monopoles and antimonopoles within small-holonomy (anti)calorons. The 
emergence of the macroscopic adjoint scalar field <\> breaks the fundamental gauge group 
SU(2) down dynamically to the subgroup U(l). Due to the adjoint Higgs mechanism, two 

(12) 

out of three gauge modes ba^ ' acquire a temperature dependent mass, while the third 
remains massless, 



A 3 

mi = m 2 = 2e(T)\(j)\ = 2e(T)\l —2- and m 3 = 0. (2.16) 

V 2irT 

Evaluating the Polyakov loop in a different (unitary) gauge gives rise to the conclusion 
that the ground state is two-fold degenerated with respect to the (broken) global elec- 
tric Zi symmetry. Thus, the electric phase is deconfining. The temperature evolution 
of the effective gauge coupling e is derived from the demand for thermodynamical self- 
consistency, it reaches a plateau value rapidly as temperature increases (T ^> T C; e) and 
diverges logarithmically for T \ T Ci e, 

e(T)cx-log(r-r CiS ). (2.17) 

(1 2) 

Therefore, the massive gauge <5a/i ' modes become infinitely heavy and decouple at T C} e- 
The ground state in the deconfining phase is composed of interacting calorons and 
anticalorons of topological charge-modulus one and trivial holonomy. Screened magnetic 
BPS saturated monopoles are spatially isolated defects in the electric phase. Screening oc- 
curs due to short-lived magnetic dipoles provided by intermediary small-holonomy LLKvB 
(anti)calorons and due to all other stable and screened (anti)monopoles. 
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Figure 2.1: Temperature evolution of the effective gauge couplings e and g as a function 
of the dimensionless temperature = The figure is taken from [23] . 

2.3 The preconfining phase 

At T Cj e = 13.87-n|s the effective gauge coupling e diverges. Thereby, magnetic monopoles 
and antimonopoles, which are generated by the dissociation of large-holonomy calorons, 
become massless and condense pairwise, thus terminating the deconfining phase. Note 
that for T < T c ^e, the average caloron-anticaloron holonomy gradually increases with 
decreasing temperature. After a spatial coarse-graining, the thermal ground state of the 
Bose condensate of interacting monopoles and antimonopoles is entirely described by a 
macroscopic complex scalar field ip and a pure gauge a° ,9S ; only gauge fields transforming 
under U(l) survive the electric-magnetic phase transition. The macroscopic complex scalar 
field (f turns out to be quantum mechanically and statistically inert. Interactions between 
monopoles mediated by pure gauges generate isolated but unstable defects. These defects 
are closed magnetic flux lines that are composed of magnetic monopoles moving oppositely 
directed to each other in the vortex core along the flux lines. The closed flux lines collapse 
as soon as they are created, thereby inducing a negative pressure. It should be noticed that 
the magnetic flux lines need to be closed due to the absence of isolated magnetic charges in 
the monopole condensate. The spatially homogeneous and BPS saturated complex scalar 
field (p breaks the dual gauge symmetry U(1) D dynamically: the stable and propagating 
excitations in the magnetic phase are massive dual gauge modes. 

In [10] equations for phase and modulus of the macroscopic complex scalar field are 
derived. The modulus of the field ip is found to be 




(2.18) 



The effective action for (p reads 




(2.19) 
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where interactions between (screened) monopoles are absent and is an externally pro- 
vided Yang-Mills scale. Interactions are accounted for in analogy to Sec. 12.21 the topolog- 
ically trivial sector a° is decomposed into 

< = < 9S + 5a° (2.20) 
and is minimally coupled to <p. The unique effective action including interaction reads 

s = f dT / A (i F - F - + + if) • (2 - 21) 

where the Abelian field strength of the dual gauge field is given by 

F^ = {d^-d v a D ^ (2.22) 
and the covariant derivative involving the effective magnetic coupling g by 

V» = dp + \ga°. (2.23) 

A pure-gauge solution a® ,3S to the equations of motions for the dual gauge-field in the 
background of <p is found. 

The evaluation of the Polyakov loop suggests that the electric Z2 degeneracy, as oc- 
curred in the electric phase, no longer exists in the magnetic phase: the ground state of 
the magnetic phase is unique and confines fundamental, heavy and fermionic test charges. 
Nevertheless, massive gauge modes still propagate because the Polyakov loop does not van- 
ish entirely. Therefore, the magnetic phase is called preconfining. The dual gauge group 
U(1)d is dynamically broken due to the emergence of the macroscopic scalar field (p. As a 
consequence, the dual gauge excitation 5a® becomes Meifiner massive via the dual Abelian 
Higgs mechanism, 

m D = g(T)<p. (2.24) 

The evolution of the temperature dependent effective gauge coupling g is predicted by 
thermodynamical self-consistency. The coupling vanishes for T /* T C) e and diverges loga- 
rithmically for T \ T Cj m- 

<?oc -log(T-T C)M ), (2.25) 

where T Cj m denotes the temperature where the transition to the center phase takes place. 
The typical energy of a non-selfintersecting center- vortex loop (CVL) is oc g^ 1 . 

During the electric-magnetic phase transition, the number of polarizations of the 'pho- 
ton' jumps from two to three, thereby inducing a discontinuity in the energy density. The 
negative pressure of the ground state arises due to an equilibrium between vortex-loop 
creation by dissociation of large-holonomy calorons and the annihilation of vortex loops by 
contraction. The non-vanishing pressure Pf^ and energy density p M of the ground state 
evaluate as 

Pm = -Pm = -^A 3 M T. (2.26) 

Across the electric-magnetic phase transition at T Cj e, where e = 00 and g = 0, the pressure 
is continuous (see Fig. 12.21) relating the scales A# and Am: 

1\ 1/3 



A E = [ - ) A M . (2.27) 
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Figure 2.2: Ratio of the energy density p and temperature T 4 across the electric-magnetic 
phase transition as a function of the dimensionless temperature Xe = The dashed 

line represents the continuation of the energy pe of the electric phase (solid black line) 
for decreasing temperature T < T c ^e (supercooled state, m D = 0). The solid grey line 
represents the energy density pu in the magnetic phase for increasing temperature (m D > 
0). As long as no additional energy is available, the system remains in a supercooled state 
until a temperature Xe = 12.15 is reached. The figure is taken from |23j. 



The magnetic phase is not detected by finite-size lattice simulations, since the mono- 
poles condensate posses infinite correlation length (oc (M m+a ) _1 ), where M m+a is the sum 
of the monopole and antimonopole mass after screening: 

87T 

M m+a = — . (2.28) 

2.4 The confining phase 

First we provide some facts on the Abrikosov-Nielsen-Olesen (ANO) vortex. When embed- 
ded in three space dimensions, a point-like two-dimensional ANO vortex becomes a vortex 
line. A mesoscopic description of a static ANO vortex is given by the action of Eq. (|2.21l) 
where the potential is absent. A BPS saturated solution to the equation of motions, fol- 
lowing from action (|2.21|) . can be found that carries one unit of magnetic flux (^) and has 
vanishing core size. Outside the vortex core the pressure P v (r), which is isotropic in the 
x%-X2 plane, reads 

^(r) = -^V, (r>0), (2.29) 



where r is the radial vector in the x\-x-i plane. Notice the minus sign on the right hand side 
of Eq. (I2.29p . For a finite energy, the length of the ANO vortex line must be finite. The 
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Figure 2.3: The oppositely directed center fluxes in the core of the intersection of a selfin- 
tersecting center-vortex loop generate an eddy where an isolated magnetic Z2 monopole is 
located. 

configuration is static as long as it possesses cylindric symmetry, but as soon as the vortex 
is bend the configuration becomes unstable: the pressure inside the vortex loop is more 
negative than outside. Thus, the vortex collapses as soon as it is created at finite coupling 
g. Notice that in the limit where g diverges the pressure vanishes. This implies that the 
formerly unstable vortex loop becomes a stable and massless particle-like excitation for 
temperatures below T Cj m- The typical core size d and energy E v of a CVL are given by 



By the collective dissociation of large-holonomy calorons and anticalorons in the pre- 
confining phase, isolated and closed magnetic flux lines start to form. At T Cj m, where the 
magnetic coupling g diverges logarithmically, the dual gauge field becomes infinitely heavy. 
Thus, a complete decoupling of the dual gauge modes takes place at the magnetic-center 
phase transition. As a result, only contact interactions between center- vortex-loops are 
possible. 

The decay of the monopole-antimonopole condensate and the subsequent formation of 
the (Bose) condensate of CVLs is described by a macroscopic complex scalar field $ in a 
potential Vc(<&). The expectation of $ is proportional to the expectation of the 't Hooft 
loop operator which is a dual order parameter for confinement. CVLs in the magnetic 
phase are created by phase jumps of $ and an increase in the modulus of <3>. This process 
continues until relaxes to one of the Z2 degenerated, energy and pressure free minima of 
a potential Vc- The phase of $ is given by a line integral of the dual gauge field a° along 
a spatial circle of infinite radius 5^ = °° measuring the (quantized) magnetic flux through 
the minimal surface Ai s n=oo. The creation of a CVL now proceeds by an infinitely thin 
flux line and its flux reversed partner traveling in from infinity and intersecting with the 
S^ 00 , thereby piercing the surface M. s r=co . The energy needed to create a single center- 
vortex loop is provided by the potential Vc- Selfintersecting and therefore massive CVLs 
come into existence when generated single CVLs that move fast enough to convert some 
of their kinetic energy into mass collide and merge, thus creating selfintersections. Each 
intersection point carries one unit of magnetic charge, see Fig. 12.31 where each sign is 
equally likely. The spectrum of excitations is equidistant since the mass of a soliton with 




and 




(2.30) 
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N self-interactions is given by NAq, Ac being the Yang-Mills scale. Modulo the charge 
multiplicities the number of distinct topologies of iV-fold selfintersecting solitons is given 
by the number of distinct topologies of connected bubble diagrams with N vertices in a 
scalar Xcp 4 quantum field theory. In Fig. 11.11 the topologies of CVLs with up to N = 3 
selfintersections are shown. If subjected to mixing with theories possessing propagating 
photons the only stable excitations are non-selfintersecting and one-fold selfintersecting 
CVLs. This is due to the repulsive or attractive forces between the charges of CVLs with 
more than one selfintersection. CVLs without selfintersection, however, are unstable in the 
presence of a noisy environment, see Sec. 12.61 

In [10] the potential Vc for the macroscopic field $ is found to be 



^=l'#-Ac*V#-Ac*l. (2.31) 



satisfying the following properties (see also |24j): 

(i) Vc is invariant under center jumps <3? — ► exp(i7r)$ only; 

(ii) it allows for the creation of spin- 1/2 fermions by a forward- and backward tunneling 
which corresponds to local center jumps of $'s phase; 

(iii) the degenerated minima of Vc have zero energy density and are related by local center 
jumps; 

(iv) a mass scale occurs to parameterize the potential Vc; 

(v) Vc needs to be real. 

The process of relaxation of $ to one of the minima of Vc is described by the action 

S = J dx 4 Hfi&d^ - \v c ^j . (2.32) 

Once $ has reached Vc's minima, quantum fluctuation <5<3? are absent because every po- 
tential fluctuation would be harder than the maximal resolution. 

Neglecting contact interactions between and internal degrees of freedom within solitons 
as well as long-range interactions between charges mediated by photons, the naive series 
for the total pressure Pc at temperature T represents an asymptotic expansion in powers 
of a suitably defined coupling coupling constant A = exp(— Ac/T). That is, the sum 

oo 

over partial pressures Pc,N of spin-1/2 states arising from solitons with N selfintersections 
seems to converge up to a critical, temperature-dependent value N C (T), but converges 
when including higher contributions. This signals that the assumption that solitons with 
arbitrary N are stable breaks down to hold for N > N c as a consequence of contact 
interactions which increase due to the higher density of intersection points and vortex 
lines. Though formally divergent, the sum over partial pressures Pc,N turns out to be 
Borel summable for negative (unphysical) values of A. The inverse Borel transformation 
is meromorphicH in the entire A-plane except for a branch cut along the positive-real axis. 
Continuation to the physical region A > leads to a sign-indefinite imaginary part which 

3 A meromorphic function is holomorphic on an open subset of the complex plane except for a set of 
isolated poles. 
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is smaller than the real part for sufficiently small temperatures. Complex admixtures to 
the pressure become manifest as turbulence-like phenomena in the plasma and thus violate 
thermal equilibrium. At zero temperature, the pressure of the ground state is precisely nil. 
Because of the over-exponential rise of spin-1/2 fermion states with increasing temperature, 
the imaginary part starts to dominate the pressure and the thermodynamical description 
of the system begins to fail (violation of spatial homogeneity). That is, at temperature 
Tfj ~ Ac, the entropy wins over the Boltzmann suppression in energy and the partition 
function diverges. This is an indication for the Hagedorn transition to the preconfining 
phase. For details see [25]. Similar behavior is observed for the expansion of the energy 
density [26] 

oo 

PC = PC,N, (2-34) 

where pc,N is the energy density of soliton states with N selfmtersections and mass NAc- 
Demanding for continuity of the negative pressure across the magnetic-center phase 
transitions yields a relation between Am and Ac>: 

A M oc 2 1 / 3 A c . (2.35) 

The question may arise whether there are stable sefhntersecting vortex-loops in the 
magnetic phase. By the decay of the macroscopic ground state in the magnetic phase its 
energy density is used to create selfmtersecting CVLs. An A^-fold twisted CVL possesses 
a mass NAc, where the Yang-Mills scale is about Ac oc Th- For this reason, the potential 
in the magnetic phase cannot provide enough energy density to create a selfintersection in 
the magnetic phase. 



2.5 The postulate SU(2) cmb = U(l) r 

We have mentioned in Sec. 12.21 that the spatial coarse-graining over the topologically non- 
trivial sector leads to the emergence of a macroscopic adjoint Higgs field which breaks 
the fundamental gauge group SU(2) down dynamically to the subgroup U(l). Thereby, 
two out of three gauge bosons become massive. At T C) e = 13.87^, where the electric- 
magnetic phase transition takes place, the mass of these two gauge field diverges and the 
massless mode remains exactly massless because radiative corrections are absent due to 
the decoupling from its heavy partners. 

Now consider the U(l)y factor of the electroweak gauge group S\J(2) W x U(l) y of the 
present Standard Model of Particle Physics (SM). In Quantum Electrodynamics (QED), 
the photon is observed to be unscreened and practically massless (m 7 < 10~ 14 eV) [27J. 
It is described by the gauge group U(l) the progenitor of which is the U(l)y factor. As 
stated above, there is only a single point in the phase diagram of SU(2) Yang-Mills ther- 
modynamics that exhibits a precisely massless gauge mode: the deconfining-preconfining 
phase transition at T Cj e- Therefore, in [TUl [23J [28J [29] , the postulate was pushed forward 
that the U(l)y factor of the electroweak gauge group is the unbroken subgroup of an SU(2) 
Yang-Mills theory with a scale comparable to the temperature of the cosmic microwave 
background (CMB) T CMB = 2.728 eV. This group is denoted SU(2) CMB . The photon 7 in 
the SM has to be identified with the massless gauge mode of SU(2) CMB . In analogy to the 



2.6. MOTION BY CURVATURE 



17 



gauge bosons of the SM, the remaining two infinitely massive and thus undetectable 
gauge modes of SU(2) CMB are denoted (m v ± = 2ecj) with e = oo at T C: e)- Furthermore, 
the average temperature of the universe T CMB is identified with the critical temperature 
T Ct E of SU(2) CMB . This fixes the only free parameter of the theory, the Yang-Mills scale A#, 
to Ae = j^jT cmb = 1.065 x 10~ 4 eV. For temperatures much above T c e, the effects of 
V ± are completely negligible, whereas for temperatures a few times of T CMB , these lead to 
a visible modification of the black-body spectrum at low frequencies (spectral gap) [29, 30J. 
The spectral gap could also provide for an explanation why old (estimated age ~ 50 million 
years), cold (mean brightness temperature ~ 20 K) and dilute (number density ~ 1.5 cm 3 ) 
clouds in between the spiral arms of the outer galaxy are composed of atomic instead of 
molecular hydrogen, and why these clouds are stable [31J. 

Regarding the transition towards the preconfining phase, the postulate SU(2) CMB = 
U(l)y implies that the photon will acquire a Meifiner mass because of the coupling to 
the newly emerging superconducting ground state (condensate of magnetic monopoles). 
The system, however, remains in a supercooled state down to T = 12.15^- due to the 
shift in energy density at T c ^e (additional degree of freedom), see also Fig. 12.21 In (23] 
an upper bound for the time the photon remains massless was estimated to be ~ 2.2 
billion years. The observed intergalactic magnetic fields can possibly be explained by 
the electric-magnetic phase transition. Conventional superconductors consist of a Cooper- 
pair condensate of electric charges and expel magnetic fields from their interior (Meifiner 
effect). If the occurrence of intergalactic magnetic fields is addressed to the emergence of 
a superconducting ground state, this leads to the conclusion that a magnetically charged 
object in the gauge group SU(2) CMB is interpreted as an electrically charged object with 
respect to U(l)y. Therefore, the ground state of the magnetic phase is a condensate 
of electrically charged monopoles and antimonopoles with respect to U(l)y, and thus 
generates intergalactic magnetic fields. 



2.6 Motion by curvature 

Here, we would like to illustrate how the curve shrinking process is induced by the curvature 
of a CVL. Recall that the vortex loop is generated by the bending of a straight ANO vortex 
line which exhibits isotropic pressure perpendicular to its symmetry axis. Now consider 
a situation where a CVL of an isolated SU(2) Yang-Mills theory is (locally) embedded 
into a flat two-dimensional surface at m D < oo and d > 0. Then, a hypothetical observer 
measuring a positive (negative) curvature of a segment of the vortex line experiences more 
(less) negative pressure in the intermediate vicinity of this curve segment (see Sec. I2.4|) 
leading to its motion towards (away from) the observer, see Fig. 12.41 The (inward directed) 
velocity of a point in the vortex core will be a monotonic function of the curvature at 
this point. On average, this shrinks the CVL. Alternatively, one may globally consider 
the limit m D — > oo, d — > 0, that is the confining phase, but now taking into account the 
effects of an environment that locally relaxes this limit (by collisions) and thus also induces 
curve shrinking. This situation is described by a curve shrinking flow in the dimensionless 
parameter r 



<9 T f(£,r) = -d|f(£,-r), 



(2.36) 
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Figure 2.4: Highly space-resolved snapshot of a segment of a center-vortex loop. The 
pressure Pi in the region pointed to by the normal vector n is more negative than the 
pressure P e thus leading to a motion of the segment along n. 



where x is a point on the planar CVL, £ is arc length, and a a string tension effectively 
expressing the distortions induced by the (noisy) environment. After a rescaling to dimen- 
sionless variables, 

x = ^fox and s = (2.37) 
flow equation (|2.36|) assumes the form: 

9 T x(s,r) = a s 2 x(s,r). (2.38) 
In the following sections, we will resort to the dimensionless flow equation. 



Chapter 3 



Mathematical Prerequisites: Curve 
shrinking flow 

In the 1970s, William Thurston developed a program for the classification of three-dimen- 
sional manifolds. It had a great impact in the field of three-dimensional topology and 
revealed a very strong connection between low-dimensional topology and differential ge- 
ometry, especially between hyperbolic geometry and Kleinian groups [32, 33J. 

Now consider a smooth closed (that is, compact and without boundary) manifold M. 
equipped with a smooth time-dependent Riemannian metric g(r). A (topological) manifold 
is a topological space which is locally homeomorphic to a Euclidean space, but with an 
generally more complicated global structure. A manifold equipped with a Riemannian 
metric g is a real differentiable manifold M., in which each tangent space is endowed with 
an inner product g in a manner that varies smoothly from point to point. It should be 
noted that not every manifold admits a geometry. The Ricci flow is a means of processing 
the metric g by the evolution of g under the following partial differential equation (PDE) 

^g(r) = -2Ric( 5 ), (3.1) 

where Ric is the Ricci curvature. In local coordinates the coefficients Rij of the Ricci 
curvature tensor are given by a contraction of the Riemannian curvature tensor R l jkh 
Rij = R k ikj- Roughly speaking, the Ricci-flow contracts regions of positive curvature and 
expands those of negative curvature, thereby smoothing out irregularities in the metric. 
In this spirit, it is formally analogous to the diffusion of heat that describes how an irreg- 
ular temperature distribution in a given region tends to become more homogeneous over 
time. An example of its application is the proof of the two-dimensional uniformization 
theorem, which states that any surface admits a Riemannian metric of constant Gaussian 
curvature. Here, the (suitably renormalized) Ricci flow is used to conformally deform a 
two-dimensional metric on Ai into one of constant curvature |34j . Richard Hamilton intro- 
duced the Ricci flow with the intention to gain insight into the geometrization conjecture 
proposed by William Thurston in 1980 [35]. The geometrization conjecture is the analogue 
for three-manifolds of the uniformization theorem for surfaces and implies several other 
conjectures, such as Thurston's elliptization conjecture or the Poincare conjecture. Let us 
first consider the Poincare conjecture, which was originally posed as a question at the end of 
an article by Henri Poincare in 1904. In its standard form, it states that every simply con- 
nected, compact three-manifold without boundary is homeomorphic to the three-sphere. 
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A more precise phrasing is that the fundamental group of a closed three-manifold M. is 
trivial, if and only if Ai is homeomorphic to the three-sphere. Now, the geometrization 
conjecture concerns the topological classification of three-dimensional smooth manifolds. 
The original phrasing of Thurston goes as follows [35]: "The interior of every compact 
3-manifold has a canonical decomposition into pieces which have geometric structures'0. 
In three dimension there are precisely eight geometric structures called the eight Thurston 
(model) geometries (involving the spherical geometry 5 3 , the Euclidean geometry R 3 , the 
hyperbolic geometry H 3 , the geometry of S 2 x R, the geometry of H 2 x R, the geometry 
of the universal cover of ST^R), the nil geometry and finally the sol geometry). The 
canonical decomposition is carried out in two steps. In the first stage, also referred to as 
the prime decomposition, one cuts a three-manifold M. along two-spheres embedded in 
M. such that neither of the obtained manifolds is a three-ball, then one glues three-balls 
to the resulting boundary components. This decomposition is unique up to the sequence 
and additional three-balls. The second stage involves cutting along certain tori that are 
nontrivially embedded in M. obtaining a three-manifold the boundary of which consists of 
tori. Hamilton's basic idea was to place an arbitrary metric g on a given smooth manifold 
M. and to dynamically deform M. by the Ricci flow to yield one of Thurston's geomet- 
ric structures. Hamilton succeeded in proving that a closed three-dimensional manifold, 
which carries a metric of positive Ricci curvature, is a spherical space form that acts like 
an attractor under the Ricci flow [36]. This is known as the Hamilton theorem. How- 
ever, in general, the Ricci flow can be expected to develop a singularity in finite time. 
Then, in a series of eprints starting in 2002, Grigori Perelman sketched a proof for the 
geometrization conjecture [37j . Thereby, Perelman modified Hamilton's program to prove 
Thurston's geometrization conjecture by stopping the Ricci flow once a singularity has 
been formed, then carefully performing 'surgery' on the evolved manifold, systematically 
excising singular regions before continuing the flow. This is called Ricci flow with surgery. 

The results obtained in this thesis heavily depend on the important work on the curve 
shortening flow done by Gage and Hamilton [38], and Grayson [MJHOj. The curve shorten- 
ing flow, also known as heat equation on immersion^], is the one-dimensional analogue to 
the Ricci flow and originally inspired Hamilton in the development of the Ricci flow. Let 
us now consider the properties of curve shrinking flows in two and three space dimensions. 



3.1 Embedded curves without selfintersection 
3.1.1 Planar curves 

Consider a family of smooth, closed curves x(s, r) of length L embedded^] in a two-dimen- 
sional flat plane R 2 , where x is a point along the curve, s G [0, L] is the arc length that 
is unique only up to a constant and r € [0, T] the flow parameter which parametrizes the 
family. The initial curve x(s,0) evolves as a function of 'time' r to x(s,r). The Euclidean 

*A geometric structure is defined to be a space modeled on a homogeneous space (X, G), where X is 
a manifold and G is a group of diffeomorphisms of X such that the stabilizer of any point x G X is a 
compact subgroup of G [35] • For every x in X, the stabilizer subgroup of x (also called the isotropy group 
or little group) is defined as the set of all elements in G that fix x: G x = {g S G\g ■ x — x} 

2 An immersion is a local embedding. 

3 An embedding is a map / : X — » Y between differentiable manifolds X and Y where the map / is a 
homeomorphism between X and its image f{X). 
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Figure 3.1: The Euclidean curve shortening flow. The arrows point towards the unit normal 
n and the length of the arrows is proportional to the curvature k. 

curve shortening flow is defined as 

<9 T x(s, r) = dgx(s, r) = k(s, r) n(s, r), (3.2) 

where the derivative d T := Jpp is taken along a fixed value of s. This is a parabolic, nonlinear 
second-order partial differential equation, where n is the inward-pointing Euclidean unit 
normal and k the scalar curvature, defined as 

k(s,r) = \d*x(s,T)\ =det(a s x(s,r),afx(s,T)), (3.3) 



with |v| = y/v ■ v, v • w denoting the Euclidean scalar product, and det(-,-) denotes the 
determinant of the 2x2 matrix created by two 2x1 vectors. It is a standard result 
for parabolic equations that solutions exist for a short time and are unique. In the curve 
shortening flow, the curve x(s, r) is deformed along its unit normal n(s, r) at a rate that 
is proportional to its curvature k(s,r). This flow deserves the attribute curve shortening, 
because its flow lines in the space of closed curves are tangent to the gradient for the curve 
length functional, see Eq. (I3.17p . For the remainder of this section, we assume that a 
solution to Eq. (|3.2|) exists on the maximal time interval [0,T). A more visual description 
of this flow is the evolution of an elastic band in a viscous medium. If the tension in 
the elastic is kept constant then its behavior is approximately determined by Eq. (|3.2p . 
see also Fig. 13.11 Since motion normal to the curve affects arc length, s is not preserved 
under curve shrinking. Thus, s and r are not independent and commute according to the 
following rule 

d T d s = d s d T + k 2 d s . (3.4) 

Therefore, we introduce the curve parameter u (modulo 2tt) related to s by 

ds = \d u x\ du. (3.5) 



The quantity \d u x\ can also be thought of as an arc length density. Actually, Eq. (|3.5|) 
defines s. The operator d s then writes as 

9 s = rJ-^du. (3.6) 
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In the following, we resort to a slight abuse of notation by using the same symbol x for 
the functional dependence on u or s. Let us now introduce coordinates in M 2 , x(u, r) = 
(x(u, t), y(u, t)) (where T denotes the transpose). The tangent vector to the curve is 
given by <9 u x, and thus we define the unit tangent t vector as 



t(u,r) := 



dux _ 1 / d u x 
\d u x\ |a u x| V d uV 

The unit normal is then given by 



n(u,r) := 



1 



-d u y 



\d u x\ V 9 u x 

The unit tangent and normal vectors are written in terms of arc length s as 



t(s,r) 



d s x 
d s y 



and 



n(s,r) 



9 s x 



(3.7) 



(3.8) 



(3.9) 



So we can write the Prenet-Serret formulas, which describe the kinematic properties of a 
point (particle) that moves along the planar curve x as 







k 
-k 



ds \ n 

where the curvature, when expressed in coordinates, is 

k( s , t) = d s x d 2 s y - d 2 s x d s y. 
The circumference of the curve L at time r is defined as 

rL{r) r2n 

L(t)= \ ds= du |d u x(u,r)| . 
Jo Jo 

The evolution of L under the flow is given by 



For the area A enclosed by the curve we have 



du \ d u x\ k 2 



Mr) ^ - 2 



rL(r) 

/ ds x(s, r) • n(s, t) 
Jo 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



Surprisingly, the time derivative of the enclosed area remains constant under curve shrink- 
ing, 



A(t) :- 



dA(r) 
dr 



-2tt. 



(3.15) 



For planar curves, the decreasing integral ds \ k\ measures the total change in angle. In 
the special case of convex planar curves, ds\k\ = J Q L ds k measures the winding number 
of the curve and is an invariant of the flow (until a singularity develops) . 
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In [4T], Grayson stated that under the flow Eq. (|3.2p 'the curve is shrinking as fast as 
it can using only local information'. Let us see how this statement can be understood. 
Consider the curve length L{r) = J^ 77 du\d u x\. To take the time derivative of L we 
differentiate |<9 u x| 2 with respect to r and obtain 

3 r |<9 u x| = — — -<9 u x • d T d u x. (3.16) 

Substituting this into L(t) and integrating by parts, we obtain the following expression 
for the rate of decrease of curve length 

f L (r) 

L(t) = - dskn-d T x. (3.17) 



J o 

Therefore, Eq. (|3.2p expresses the local condition that the rate of decrease of L{t) is 
maximal with respect to a variation of the direction of the velocity <9 T x of a given point 
on the curve at fixed magnitude |<9 r x| [12]. However, the magnitude |<9 r x| is not in general 
the speed which maximizes L(t). 

Setting A(t = 0) = Aq, the solution to Eq. (I3.15P is 

A(t) = A - 2vrr. (3.18) 

By virtue of Eq. (|3.18|) the critical value T, where A and with it the curve vanishes, is 
related to Aq as 

T=^. (3.19) 

The isoperimetric ratio is defined as and the isoperimetric inequality states that 
L{rf 



A(r) 



> 4tt. (3.20) 



Equality is achieved if and only if the curve is a circle. Therefore, one can consider it a 
measure of 'how circular' the curve is. 

In 1983, Gage showed that when a smooth convex curve evolves according to Eq. (|3.2p . 
the isoperimetric ratio decreases, so that if A — > 0, then L — > and the curve shrinks 
to a point [13]. In 1984, Gage showed that a convex curve is becoming circular and ^ 
approaches 47r, as the enclosed area approaches zero, provided that the curvature does not 
blow up prematurely, that is the curve does not form a cusp [M]. As a consequence, the 
ratio of the circumscribed ratio to the inscribed ratio converges to unity. This can be 
considered a C°-convergence to the circle. Hence, in the absence of singularities, a strictly 
convex and embedded curve remains convex and embedded under the evolution. 

In 1986, Gage and Hamilton showed that for convex curves the curvature does not blow 
up prematurely for lim T A[t) = (38) ■ Thus, the curve remains convex and becomes 
circular, as it shrinks to a point for r f T, where < T < oo. The curve shrinks to a 
circle in the sense that: 

(i) the ratio approaches unity; 

(ii) the ratio of the maximum curvature to the minimum curvature x 21 ^ approaches unity 

^min 

(C 2 -convergence) ; 
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(iii) the higher order derivatives of the curvature k converge to zero uniformly (C°°- 
convergence) . 

In 1987, Grayson showed that embedded (non-selfintersecting) planar curves become con- 
vex before T without developing singularities [H]. Thus, this completes the proof of the 
well known Gage-Hamilton-Grayson theorem that curve shortening determined by Eq. (|3.2p 
shrinks embedded plane curves smoothly to points, with round limiting shape. It is impor- 
tant to note that some planar curves, which are immersed but not embedded will surely 
develop singularities, e.g. the figure-eight of Sec. 13.21 or the Limagon of Pascal. 

Consider the set of all Euclidean transformations in M 2 , that is the set of all rotations, 
translations and reflections of a figure in M 2 . Such a transformation ET : M. 2 — > M 2 is a 
function of the form 



where U is an orthogonal 2x2 matrix and a € R 2 . The Euclidean curve shortening flow 
is defined in terms of the Euclidean curvature k and the Euclidean unit normal n that are 
invariant under Euclidean transformations ET. 

3.1.2 Space curves 

There are several possibilities of generalizing the curve shortening flow. One is the mean 
curvature flow, which is the generalization of Eq. (|3.2|) for hyper surf aces. In this case, the 
results of Sec. 13. l.fl continue to hold for convex curves, but for non-convex curves they do 
not [15]. 

For our purposes, it is more interesting to look at the extension of the curve shortening 
flow for curves embedded in the three-dimensional Euclidean space M 3 . Consider a contin- 
uous, differentiable (and not necessarily closed) space curve x(s,r) embedded in M 3 . The 
tangent, normal and binormal unit vectors are denoted t, n and b, respectively, also called 
Frenet-Serret frame, and defined as follows: 

t is the unit vector tangent to the curve, pointing in the direction of motion: t = d s x; 
n is the normalized derivative of t with respect to the arc length s of the curve: n = jfjfr; 
b is the cross product of t and n: b = t x n. 

The Frenet-Serret formulas for a point on the space curve are given by 



where k is the curvature and t the torsion. The Frenet-Serret formulas effectively define 
the curvature and torsion of a space curve. It should be noted that the existence of a 
Frenet-Serret frame requires |fcn| 2 > 0. That is, a particle traveling along the curve must 
experience acceleration. The evolution equation for space curve assumes the same form as 



ET(x) = Ux + a, 



(3.21) 




(3.22) 



Eq. 



9 T x(s,r) = k(s, r)n(s, r). 



(3.23) 



The unit normal n is not always defined, though kn always makes sense. It was shown 
by Altschuler and Grayson that solutions to the space curve flow exist until the curvature 
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becomes unbounded. However, space curves may not remain embedded in general, and sin- 
gularities will develop in the case of closed curves. A phenomenon of space curve evolution 
is that inflection points (k = 0) may develop during a time interval on which the curvature 
is bounded. When this happens the curvature becomes zero and the torsion infinite at 
a point. Nevertheless the curve remains embedded since the flow ignores these types of 
singularities in the torsion [46JII7]. A rather surprising property of space curve evolution 
is, that the formation of a singularity is a planar phenomenon. A space curve is said to 
be planar at a point (s',t') if the ratio of torsion and curvature vanishes, t(s',t') = 0. 
In [46], Altschuler showed that if a space curve develops a singularity at (s',r'), then 
]imfg T \_+r g t )T r\ f(s,r) = 0. Furthermore, Altschuler showed that the space curve is either 
asymptotic (r — ► oo) to a planar solution which moves by homothety (self-similarity), or 
a rescaling of the solution along the singularity converges in C°° to a limiting solution 
x(s,r = oo) [46], where x(s,r = oo) is the family of planar, convex curves. The most 
trivial case of a curve moving by homothety is the circle shrinking down to a point. It 
should be noted that the conjecture due to Grayson, that singularity formation is a planar 
phenomenon, can be proven without using the language of rescalings. 



3.2 Immersed curves with one selfintersection 

When a closed curve immersed in a plane evolves by its curvature according to Eq. (|3.2p . it 
remains smooth until its curvature blows up. From Sec. 13.1.11 we know that an embedded 
closed curve cannot develop a singularity until it shrinks to a point, where the limiting 
shape of the curve converges in C°° to circle. In marked contrast to this behavior, it was 
shown by Grayson that an immersed curve can evolve by the curvature flow such that its 
area vanishes, but its isoperimetric ratio converges to oo. Such a curve, namely a figure- 
eight, was investigated in [40]. A figure-eight is the simplest non-embedded curve and is 
defined to be a smooth immersion into the plane with exactly one double point, and a total 
rotation number zero, 

L 

dsk = Q. (3.24) 

Here, s is arc length, L the curve length and k the scalar curvature. Such a curve divides 
the plane into three disjoint areas two of which are finite and denoted (the unsigned areas) 
A\ and A%. Let x(s,r = 0) be figure-eight which evolves to x(s,t) according Eq. A3. 2D 
for < r < T. The curvature is unbounded as r — > T. The main result of [30] is that 
the isoperimetric ratio ^- converges to oo as r — > T if and only if the loops bound regions 
of equal area, Ai(0) = ^(0). This in turn implies that ^- for a curve with unequal-area 
loops is bounded as r — » T. 

Since for immersed curves the number of double points is a non-increasing function of 
time [48], a figure-eight remains a figure-eight until one of its loops collapses or the flow 
encounters a singularity. The curve stays smooth and the flow continues until A\ or A^ 
converge to zero. For the total area A of a figure-eight we have 

A(t) = A 1 (t)+A 2 (t). (3.25) 

The time derivative of the area enclosed by one of loop of the curve is equal to — | J ds k\ 
over the loop. Unlike the case of a non-selfintersecting curve, the rate of change of the 
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total area is not longer constant, but constrained as 



-4tt < 



- dr ~ 



< -2tt. 



(3.26) 



However, we have the nice property of figure-eights that the difference of areas bounded 
by the two loops of x(s,r) remains constant under the flow evolution: 



Aside from a number of applications in differential geometry, curve shortening flows 
are also used in multi-agent systems, such as mobile autonomous robots [12], in image 
processing where the flow provides an efficient way to smooth curves representing the 
contours of objects, or in computer vision. For a complete account of many of the results 
of curve shrinking see [451 149] . 

In the following, we suppress the functional dependence on u in the argument of x and 
n and write x(r) := x(u, r) and u(u, r) := n(r) 



A\(t) — A2(t) = const. 



(3.27) 



Chapter 4 

Non-selfintersecting center- vortex 
loops 

We apply curve shrinking to the N = sector in the sense of Sec. 12.61 It should be 
noticed that the restriction of the motion of a CVL to a two-dimensional flat plane is a 
major assumption which needs to be supplemented by additional physical arguments for 
its validity. 

4.1 Wilsonian renormalization-group flow 

In this section, we exploit the concept of renormalization-group transformations to yield 
an effective 'action' that enables us to compute statistical quantities. The renormalization 
group allows one to investigate the change in the physical parameters of a system which is 
associated with the change in scale (energy or resolution) and necessary to keep the physics 
constant. In our case, the change of scale corresponds to a change of the resolution Q used 
to probe the system. Here, the resolution Q is a strictly monotonic decreasing function 
of the flow parameter r. The change in parameters of the effective 'action' is implicitly 
determined by a renormalization-group flow in r. 

4.1.1 Geometric partition function 

Let us now interpret the process of curve shrinking determined by Eq. (|3.2p as a renorm- 
alization-group transformation of a statistical ensemble made up of planar N = CVLs. 
A partition function, which is the sum over suitable defined weights of the members in the 
ensemble, is considered to be invariant under a decrease of resolution Q determined by the 
flow parameter r. Physically, r is a monotonically decreasing function of Q/Qo, where Q 
(Qo) are mass scales associated with the actual (initial) resolution applied to the system. 
The role of Q can also be played by the finite temperature of a reservoir that is coupled 
to the system. 

To define a suitable weight, we devise an ansatz for the effective 'action' S = S[x(t)] 
in geometric terms of the curves in the ensemble, since these are the only accessible quan- 
tities in the system of isolated non-interacting CVLs. The 'action' as a functional of x is 
expressible in terms of integrals over local densities in s. Furthermore, we take advantage 
of the following symmetries the action should possess: 
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(i) scaling symmetry x — > Ax, A G M + : for both conformal limits, A — > oo and A — > 0, where 
the curves at fixed L gets unobservable since AL — > oo and AL — > 0, the 'action' 5 should 
be invariant under further finite rescalings (decoupling of the fixed length scale a -1 / 2 ); 

(ii) Euclidean point symmetry in M 2 , that is the group of all rotations, translations and 
reflections of a figure (curve) in the plane: sufficient but not necessary for this is a repre- 
sentation of S in terms of integrals over scalar densities with respect to these symmetries. 
That is, the 'action' density should be expressible as a expansion in series involving prod- 
ucts of Euclidean scalar products of JprX, n £ N + , or constancy. However, scalar integrals 
can be constructed which involve non-scalar densities. For instance, consider the area A 
enclosed by curve and given by 

-L(r) 



AM = i 



/ ds x(t) • n(r} 
Jo 



(4.1) 



The density x • n in this expression is not a scalar under translations. 

We now decompose the effective 'action' into a conformal and a non-conformal factor 

S = F c x F nc , (4.2) 

where in addition to Euclidean point symmetry F c is invariant under x — ► Ax, whereas F nc 
is not. In principle, infinitely many operators can be defined to contribute to F c . Since 
the evolution generates circles for r /" T and thus homogenizes the curvature, higher 
derivatives of k with respect to s rapidly converge to zero [38]. We expect this to be true 
also for Euclidean scalar products involving higher derivatives JprX. To yield conformally 
invariant expressions such integrals need to be multiplied by powers of \[~A and/or L or the 
inverse of integrals involving lower derivatives. At this stage, we are not able to constrain 
the expansion in derivatives by additional physical or mathematical arguments. To be 
pragmatic, we simply set F c equal to the isoperimetric ratio: 

(«) 

We consider the non-conformal factor F nc in S as a formal Taylor expansion in inverse 
powers of L or A due to the conformal invariance of the curve for L, A — > oo and L, A — > 0. 
Since the renormalization-group evolution of the effective 'action' is driven by the curve 
shortening flow of each member in the ensemble, we allow for an explicit r dependence of 
the coefficient c of the lowest nontrivial power \. The idea is to include the contribution 
of higher-order operators, that do not exhibit an explicit r dependence, into a resolu- 
tion dependence of the coefficient of the lower-dimensional operators. Thus, we make the 
following ansatz 

Fnc{r) = 1 + . (4.4) 

The initial value c(r = 0) is determined from a physical boundary condition such as the 
mean length L at r = which determines the mean mass m of a = CVL as fh = aL. 
We have also considered a modified factor 

F nc (r) = 1 + ^ (4.5) 
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in the ansatz for the 'action' in Eq. (|4.2|) . 

For later use, we investigate the behavior of F nc (T) for r /" T for an ensemble consisting 
of a single curve only and require the independence of the 'partition function' under changes 
in t. Using Eq. (I3.18P in the vicinity of r = T, where the limiting of curve is a circle with 
radius R, we have 

L{t) = 2-kR = V8ir VT-t. (4.6) 

Since F c (t /T) = An, independence of the 'partition function' under the flow in r implies 
that 



c(t) (xVT^t. (4.7) 

That is, F nc approaches a constant value for r /" T which brings us back to the conformal 
limit by a finite renormalization of the conformal part F c of the effective 'action'. In this 
parametrization of S, the coefficient c(r) can thus be regarded as an order parameter for 
conformal symmetry with a mean-field critical exponent. 

4.1.2 Effective 'action' 

We now want to derive an effective 'action' 5[x(r)] resulting from a partition function Z 
for a nontrivial ensemble E. The partition function Zm is defined as the average 

M 

Z M = ^exp(-5[x i (r)]) (4.8) 

i=l 

over the ensemble E = {xi,... xm}. Em denotes an ensemble consisting of M curves 
where Em is obtained from Em-i by adding a new curve xm(«,t). We are interested in 
a situation where all curves in Em shrink to a point at the same value r = T. Because of 
T = ^4o/(2vr), we demand that at r = all curves in Em have the same initial area A$. The 
effective 'action' S in Eq. (|4.2|l (when associated with the ensemble Em we will denote it as 
Sm, and the corresponding coefficient cm) is determined by the function cm(t), compare 
with Eq. (|4.4p . the flow of which follows from the requirement of r-independence of Zm' 



-^Z m = . (4.9) 

This is an implicit, first-order ordinary differential equation for cm(t), which is in need for 
an initial condition cq,m = cm(t = 0). An obvious choice of initial condition is to demand 
that the statistic mean length L(r), defined as 

1 M 

Lm(t) = — -- 5>[xi(r)] exp(-S M [x l (r)]), (4.10) 



coincides with the geometric mean length Lm(t) defined as 

M 

M 



1 M 

£ M(T) = -J>[ Xi (T)] (4.11) 



i=l 
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at r = 0: 



L M (0) =L A /(0). 



(4.12) 



Prom this initial condition a value for co,m follows. In the case of the modified 'action' in 
Eq. (|4.5p . the choice of initial condition Lm{j = 0) = Lm(t = 0) leads to F nc (r) = which 
is equivalent to a uniform distribution. This is because initial condition (|4.12l) is identically 
fulfilled for the modified 'action' if c(0) = — Aq is chosen, then setting c(r) = — A(t) solves 
dZM/dt = trivially. While the geometric effective 'action' is thus profoundly different for 
such a modification of F nc (r), physical results such as the evolution of the variance of the 
position of the 'center of mass' agree remarkably well, see Sec. 14.2.41 We conclude, that 
the geometric effective 'action' itself has no physical interpretation in contrast to quantum 
field theory and conventional statistical mechanics where the action in principle is related 
to the physical properties of a given member of the ensemble. Rather, going from one 
ansatz for Sm to another describes a particular way of redistributing the weight in the 
ensemble which seems to have no significant impact on the physics. 

4.2 Results of simulation 

4.2.1 Preparation of ensemble 

For the curves depicted in Fig. 14.11 we make the convention that Aq = 2ir x 100. It 
then follows that T = 100 by virtue of Eq. (|3.19|) . Furthermore, we have prepared the 
ensembles such that the position of 'center of mass' (COM) coincides with the origin. It 
should be recalled that such a translation does not alter the effective 'action' (Euclidean 
point symmetry). Also note that we use the same notation Em for the primed and the 
unprimed ensemble. In Fig. 14.21 the evolution of two different initial curves under curve 
shrinking is shown. 

4.2.2 Numerical procedure 

The initial curves depicted in 14.11 are generated as follows. First, we chose a list of points 
in the x-y plane such that the initial curve x(it, r = 0) consecutively passes the points, 
where initial and final point coincide. A spline, which is made of piecewise third-order 
polynomials with C 1 continuity, interpolates each point in the list. Since Mathematica's 
SplineFit sets the second derivatives of the spline at the endpoints to zero, the first two 
points are appended to the end of the list. The analogue holds for the last two points of 
the original list. To yield a smooth curve with more than C 1 continuity each coordinate 
of the curve is fitted by trigonometric functions up to order eight in its Fourier-expansion. 
Area, length and centroid of the initial curve are computed numerically with NIntegrate, 
where the latter is given by 



Now we can prepare the ensembles as described in Sec. 14.2.11 

To simulate the flow evolution of the initial curves one searches for solutions to the 
second-order partial differential equation 




(4.13) 



d T x(u, t) 



8 u x(u,t)\ u |<9 m x(u,t)| 



<9 u x(li,T) 



(4.14) 
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Figure 4.1: Initial curves contributing to the ensembles Em- The positions of the 'center 
of mass' coincide with the origin, and all curves have the same area 200 it. 

subject to periodic boundary conditions in the curve parameter, x(u = 0, r) = x(u = 
2-7t,t), and for the initial conditions x(u,t = 0) depicted in Fig. 14.11 This is done using 
the Numerical Method of Lines. This is a technique for solving PDEs by discretizing 
in all but one dimension, and then integrating the semi-discrete problem as a system of 
coupled ordinary differential equations (ODEs) or differential-algebraic equations. Here, 
we partially discretize the flow equation Eq. (14. 14ft on a uniform grid in the parameter 
u yielding an ODE initial value problem in r that was solved by the ODE integrators in 
Mathematica's NDSolve. Fig. l4.2l indicates why this technique is called the method of lines. 
As one can also see from Fig. 14.21 a set of discrete points on the curve, although remaining 
equidistant in u, may evolve under the flow such that the spatial distances between adjacent 
points falls below the numerical precision. The flow then encounters a purely numerically 
and thus virtual singularity (not to be confused with the earlier mentioned non-virtual 
singularities at r = T). Therefore, the execution of NDSolve is broken up into several 
basic steps which are carried out separately. These steps are: 

(i) equation processing and method selection, 

(ii) method initialization, 

(iii) numerical solution, 

(iv) solution processing. 

The low-level functions that are used in Mathematica to break up these steps are ND- 
Solve' ProcessEquations (i,h), NDSolve' Iterate (iii) and NDSolve* ProcessSolutions (iv). 
NDSolve* ProcessEquations classifies the differential system into an initial-value prob- 
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100 




X 



10 



(a) 



(b) 



Figure 4.2: Plots of the evolution of planar N=0 CVLs (curve 2 and 6 of Fig. 14.11) under 
the curve shortening flow. The thick central lines depict the trajectories of the 'center of 
mass' (see Sec. I4.2.4|) which coincides with the origin at r = 0. The flow is started at r = 
and stopped at r = 100. 

lem, boundary-value problem, differential-algebraic problem, partial-differential problem, 
etc. It also chooses appropriate default integration methods and constructs the main 
NDSolve v StateData data structure. NDSolve v Iterate advances the numerical solution. 
The first invocation initializes the numerical integration methods. NDSolve' ProcessSo- 
lutions converts numerical data into an InterpolatingFunction to represent each solution. 
More precisely, the curve parameter range is divided into n equidistant intervals yielding 
n points on the curve which are generally not equidistant in space. For our simulation the 
number of points n is chosen between 130 and 300. The discretization of flow equation 
(I4.14p with respect to the variable u needs to convert the derivatives into finite differences. 
The second-order centered (with respect to to the set of sample points around x{ui)) for- 
mula for the first derivative is given by 



where h is the grid spacing. Here, finite differences of sixth order are used which are com- 
puted with Mathematica's NDSolve* FiniteDifferenceDerivative. In the following, every 
quantity involving derivatives evaluated on a discrete set of data points is computed using 
Mathematica's NDSolve~ FiniteDifferenceDerivative. After NDSolve* ProcessEquations is 
invoked the first time at r = t% = 0, the numerical solution is advanced using ND- 
Solve* Iterate by a unit 'time' step At = 1 up to r^. Then the computation is interrupted 
to compute an error estimate that indicates whether a virtual singularity is starting to 
evolve. The error estimate exploits that A{r) = Aq — 2ttt and is computed as 




x(u i+ i) - x(Uj-l) 
2h 



+ 0(h 2 ) 



(4.15) 



10 4 x (A(t 2 ) - A(n) + 2vr(r 2 - n)) 



(4.16) 
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where A{t) is given by the discrete version of Eq. (|3.14|) evaluated on the point grid given 
by NDSolve. 

Until T2 reaches T, the solution is advanced step by step as long as the error estimate 
does not exceed the empirically found value of 2. But if it does, the by then obtained 
solution is fitted at ti — 1 in such a way that a new discretization yields (spatially) well 
separated points to restart the procedure. In Fig. I4.2f b) such a situation is shown. The 
fitted curve is obtained as follows. At t%, one determines the minimal arc length s m i n which 
is the least of all arc length between adjacent points on the curve. Then, at T2 — 1, all 
those points on the curve are dropped the arc length of which to their next neighbors is 
less than the minimal arc length s m ; n . The remaining points are fitted by trigonometric 
functions, where the order of the fit is chosen to depend on r (since the curve is getting 
smoother with increasing r). In the case of the error estimate of the fitted curve exceeding 
the tolerance, the number of grid points has to be increased or the initial curve needs to be 
smoothed slightly. In order to avoid discontinuities in the r-evolution of L, A and x COM , 
and singularities in their derivatives that occur since the fit procedure generates piecewise 
defined functions, and since after the fit the values of A and L slightly deviate from their 
former values, these quantities are interpolated by polynomials for < r < T using Find- 
Fit. To improve the accuracy of L near the critical value T, the isoperimetric ratio hr is 

fitted instead of L, and L is calculated from iJ ~A. 

T 2 

The analytical results of Sec. 13.11 such as the convergence of to Air, the constancy 
of A or the vanishing of L and A for r f T are numerically well reproduced, thereby 
confirming the validity of the simulation. 

The implicit first-order differential equation ^ = for the coefficient is solved using 
NDSolve. If not set at will, the initial condition Co for c(r) was derived from Eq. (|4.12j) 
using Mathematica's FindRoot. The variance of the position of COM was computed. The 
square of the coefficient c(r) associated to the non-conformal factor was fitted with function 

c(r) 2 = k(T -rr, (4.17) 

where k and a are fit parameters. We have determined the critical exponent of the coeffi- 
cient to t| = 0.5 as r — ► T, in accordance with the theoretical value of Eq. |4~71 For checking 
purpose, we have also used Tq as fit parameter, yielding excellent agreement within the 
numerical precision. 

A CD-ROM containing the used Mathematica notebooks is attached to the thesi^]. 

4.2.3 Renormalization-group invariance of partition function 

The function c 2 m {t) is plotted in Fig. 14.31 According to Fig. 14.31 it seems that the larger the 
ensemble the closer c 2 m (t) is to the evolution of a single circle of initial radius R = yJ~^f- 
For growing M the function c 2 m {t) approaches the form 

) = k M (T-r), (4.18) 



„2 

'as,/ 



lr Tlie results for non-selfintersecting curves were obtained using Mathematica version 6.0.2 or below. 
Here, a remark concerning the used Mathematica version is in order. Due to incomprehensible reasons 
version 6.0.3 is not capable to solve the implicit ODE for the coefficient c(r), not even in the trivial case 
of an ensemble consisting of a single curve. In the case of one-fold selfintersecting curves, version 6.0.3 still 
works and was used. 
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Figure 4.3: The square of the coefficient cm(t) entering the effective 'action' Sm 



A(r) 



1 + 



Cm(t) 

W) 



for various ensemble sizes M 



■ 1, . . . , 12. Notice the early onset of the 
T does not depend on Cq m = c m (t = 0) 
and thus not on the initial choice of Z, but only on the specific choice of curves included 
in the ensemble. 



linear drop of c M (r 



The slope of c 2 m {t) near r 



where the slope ku depends on the strength of deviation from circles of the representa- 
tives in the ensemble Em at r = 0, that is, on the variance ALm at a given value Aq. 
Physically speaking, the value r = is associated with a certain initial resolution of the 
measuring device (the strictly monotonic function r(Q), Q being a physical scale such as 
energy or momentum transfer, expresses the characteristics of the measuring device and 
the measuring process), the value of Aq describes the strength of noise associated with the 
environment (^4o determines how fast the conformal limit of circular points is reached), 
and the values of co,m and /cm, see Eq. (|4.18|) . are associated with the conditions at which 
the to-be-coarse-grained system is prepared. Notice that this interpretation is valid for the 
'action' 



Sm 



A{r) 



1 + 



CM\J, 

L(t) 



only. 

If we relax initial condition Lm(0) = Zm(0) for co,m and set the initial value for co,m 
at will, the coefficient starts at the given value and rapidly adapts to the evolution depicted 
in Fig. 14.31 and respectively, Fig. 14.41 In Sec. 14.1.21 we have argued that for the modified 
'action' and the initial condition Lm(0) = Lm(0) the curves are uniformly distributed. 
Relaxing this initial condition in the case of the modified 'action', means that the curves 
are no longer uniformly distributed for r = 0. However, the uniform distribution is restored 
rapidly as the curves evolve under the flow. 



4.2. RESULTS OF SIMULATION 



35 



Cm 


E, 


E 2 


E 3 


E 4 






E 5 


E 6 


E 7 


E 8 


-m 


E 9 


Eio 


En 


E12 










z 



Figure 4.4: The coefficient ca/(t) entering the effective 'action' Sm = iTfr)" ^ 
for ensemble sizes M = 1, . . . , 12. 



4.2.4 Variance of mean 'center of mass' 

Having obtained the coefficient in the non-conformal factor of the effective 'action', we 
are now able to compute the flow of an 'observable', such as the COM position in a given 
ensemble and its statistical variance. The COM position x C q M of a given curve x(s, t) is 
defined as 

1 f L ^ 

Xoom(t) = (x om(t),2/cqm(t)) = -r-—- / dsx(s,r). (4.19) 

L \ r ) Jo 

We will present below results on the statistical variance of the COM position. 

At r = 0, the statistical variance in the position of the COM is prepared to be nil, 
physically corresponding to an infinite resolution applied to the system by the measuring 
device. In Fig. 14.51 the flow of the COM position corresponding to the initial curves 
depicted in Fig. 14.11 is shown. 

The mean COM position x COM over the ensemble Em is defined as 

1 M 

Xcom(t) = (x COM (T),y COM (T)) T = — V"x COMii (r)exp(-S , M[xi(r)]) . (4.20) 

i=i 

The scalar statistical deviation COM of 

x com over the ensemble Em is defined as 



Am,com(t) = \J varM,coM;x(T) + var^ooM^^) , (4-21) 

where 



M 



var M ,coM;x = 7j— X] ( x coM,i( T ) ~ x C om(t)) 2 exp (-S M [*i(T)]) 

i=l 

1 M 

= -4oM^) + ^^ x coM,i(^) exp(- < S A/ [x i (r)]), (4.22) 
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x 



Figure 4.5: Flow of the positions of the 'centers of masses' for the initial curves depicted 
in Fig. |4~T1 

and similarly for the coordinate y. In Fig. 14.61 plots of Am,com( t ) are shown when 
^m,com(7") is evaluated over the ensembles E%, . . . ,E\2 with the 'action' 



Sm 



L{jf 
A(t) 



1 + 



cm(t) 
L(t) 



and subject to the initial condition Lm(t = 0) = Lm{j = 0). In Fig. 14.71 the according 
plots of Am,com(''~) are depicted as obtained with the modified 'action' 



Sm 



L{jf 
A(r) 



1 + 



cm(t) 
A{t) 



and subject to the initial condition Luij = 0) = Lm{t = 0). In this case, one has 
cm(t) = — A(t) leading to equal weights for each curve in Em- Note that the slight 
qualitative deviation of the last graph E12 for small values of r regarding to the previous 
graphs in Fig. 14.61 is due to the fact that the curves which were added to the ensemble at 
last are the most twisted ones. Graph E12 still saturates at a finite value of r, nevertheless. 
The fluctuations in graph E\ of Fig. 14.61 are within the range of the numerical precision. 



4.2.5 Quantum mechanical versus statistical uncertainty 

In view of the results obtained in the last section, we would say that an ensemble of evolving 
planar CVLs in the N = sector qualitatively resembles the Quantum Mechanics of a free 
point particld! of mass m in one space dimension x. Namely, an initially localized square 



of the wave function if) with \tp(T = 0, x)\ 2 oc exp 



where Ax(r = 0) = ao, according 



2 It is no relevance at this point whether this particle carries spin or not. 



4.2. RESULTS OF SIMULATION 



37 



Am, com 


E, 

/ . . . . 1 - 


E 2 


E 3 

1 


E 4 


E 5 


E 6 


' E 7 * 


E 8 




" E 9 " 


Eio 


En 


"e u " 










t 



Figure 4.6: Plots of Am,com(t") for M = 1, . . . , 12 when evaluated with the 'action' Su = 
fl + °l(t) ) • Notice the rapid generation of an uncertainty in the COM position 
under the flow and its saturation when approaching the conformal limit t /* T. There 
also is a saturation of this limiting value with a growing ensemble size. 
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Figure 4.7: Plots of Am, C omW for M = 1, . . . , 12 when evaluated with the modified 'action' 
^ M = L A(t) \A C A(t) ) • This corresponds to a uniform distribution when evaluated with 
initial condition (|4.12|) . Notice the qualitative agreement with the results displayed in 
Fig. SSI 
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to unitary time evolution in quantum mechanics evolves as 
|^(r,x)| 2 = | exp 

2 

where H = is the free-particle Hamiltonian, p the spatial momentum, and a(r) = 



ao\Jl + (^^ij • In agreement with Heisenberg's uncertainty relation, one has during the 
evolution that 



AxAp = 




(4.24) 



The time evolution in a quantum mechanical system and the process of lowering the reso- 
lution in a statistical system describing planar CVLs share the same property: the 'time' 
(resolution) evolution generates out of a small initial position uncertainty (corresponding 
to a large initial resolution Ap) a larger position uncertainty as 'time' increases (resolution 
decreases). Possibly, future development will show that interference effects in Quantum 
Mechanics can be traced back to the non-local nature of the degrees of freedom (CVLs) 
entering a statistical partition function. 



Ht 



V>(t 



0, x)\ 2 oc exp 



(x 



(4.23) 



Chapter 5 

Selfintersecting center- vortex loops 



Let us now turn to the case of N = 1 CVLs. We proceed as far as possible in close analogy 
to the N = sector. 

5.1 Wilsonian renormalization-group flow 
5.1.1 Geometric partition function 

As in the N = sector, we interpret curve- shrinking as a Wilsonian renormalization-group 
flow. The partition function is now defined over an ensemble of JV = 1 CVLs, and we 
consider it to be independent under a change of resolution Q and thus independent of r. 
We express the effective 'action' in terms of integrals over local densities in s, and demand 
the following symmetries in order to conceive an ansatz for the effective 'action': 

(i) scaling symmetry x — > Ax, AG M + : for A — > oo, implying XL — > oo at fixed L, the 
'action' S should be invariant under further finite rescalings (decoupling of the fixed length 
scales cr -1 / 2 and A -1 ). 

(ii) Euclidean point symmetry of the plane: this is sufficiently satisfied for a representation 
of S in terms of integrals over scalar densities with respect to these symmetries. Thus, 
we can represent the 'action' density as a series involving products of Euclidean scalar 
products of ^ttX , n £ N + , or constancy. 

As in Sec. 14.1.11 we resort to a factorization ansatz as 



where in addition to Euclidean point symmetry F c (F nc ) is (is not) invariant under x — > 
Ax. In principle, infinitely many operators can be defined to contribute to F c . Since 
the evolution homogenizes the curvature, except for a small vicinity of the intersection 
point where one or both loops of the curve vanish, higher derivatives of k with respect 
to s should not be of importance. This should also hold for Euclidean scalar products 
involving higher derivatives ;JjttX. Conformally invariant expressions are obtained from 
such integrals if multiplied by powers of VA and/or L or the inverse of integrals involving 
lower derivatives. The conformal factor F c is set equal to the isoperimetric ratio, 



S = F c xF, 



nci 



(5.1) 



F c (r) 




(5.2) 



Mr) ' 
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The property of conformal invariance for L, A —* oo suggests to express the non-conformal 
factor F nc as a formal expansion in inverse powers of L or A = A\ + A%, We allow for 
an explicit r dependence of the coefficient c of the lowest nontrivial power j or -j. In 
principle, this sums up the contribution to .F nc of certain higher-power operators which do 
not exhibit an explicit r dependence. 

We restrict to the following two ansatze for the non-conformal factor in Eq. (|5.1|) . 



Fnc(r) = 1 + ^ . (5.3) 

and for the modified 'action' 

Kc(r) = 1 + ^ • (5.4) 

The initial value c(r = 0) is determined from the physical boundary condition such as the 
mean length L at r = 0. Although the modified ansatz (|5.4p in F nc of the geometric 'action' 
is profoundly different physical results such as the evolution of entropy or the variance of 
intersection of a given ensemble agree remarkably well, see Sec. 15.21 

5.1.2 Effective 'action' 

The effective 'action' Sm[x(t)] results from a partition function Zm which is defined as 
the average 

M 

Zm = y~] exp (-S M [xj(r)]) (5.5) 

i 

over the nontrivial ensemble Em = {xi,... xj\/}. The ensemble Em, consisting of M 
curves, is obtained from Em-i by adding a new curve xjf(r). The effective 'action' Sm 
in Eq. (|5.ip is determined by the function cm(t), the flow of which follows from the 
requirement of r-independence of the partition function: 

~t~Zm = 0. (5.6) 
dr 

As in Sec. l4.1.2l we obtain the initial condition co,m = cm(t = 0) to this implicit first-order 
ordinary differential equation by the constraint that the geometric mean coincides with the 
statistic mean at r = 0, 

L M (0)=L M (P). (5-7) 



5.2 Results of simulation 
5.2.1 Preparation of ensembles 

Similar to Sec. 14.2.11 all curves are normalized curves to have the same initial total area 
Aq = j4o,i + Aq : 2 and since we are now interested in the position of the intersection where 
the (anti)monopole is localized (see Sec. 12.41 or Fig. I6.ip . we have applied a translation 
to each curve in the ensemble Em such that the location of the intersections initially 
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Figure 5.1: Initial curves Xj(ii,r = 0) contributing to the ensemble Em=w- The intersec- 
tion points x- mtj i(r = 0) coincide with the origin, and all curves have the same area 200 tt. 
By definition Em=w is T-ordered. 

coincide with the origin. Again, such a transition does not alter the effective 'action' due 
to Euclidean point symmetry. 

We order the members of the maximal-size ensemble Em=w into sub-ensembles Em <i6 
such that Tj = i > T = 2 > ••• > Tm, because the critical value T of the flow parameter 
r varies from curve to curve. These ensembles Em are referred to as T-ordered. We 
have also performed all simulations with ensembles E' M<ie the members of which are 
picked randomly from Em=w and have obtained similar results for ensemble averages of 
'observables' using Em<w and E'm<ig f° r the r evolution to the left of r = min{T|xj G 
E' M<16 }. The main difference is that the computation of the coefficient, and with it the flow 
of 'observables', terminates at a smaller r since the ensembles E' M are no longer T-ordered. 

The maximal-size ensemble Em=w at r = is depicted in Fig. 15.11 with the universal 
choice Aq = 200 tt. The curves in Fig. 15.11 are arranged in a T-ordered way. We have 
Tj = i = 65 > T2 > • • • > Tm = 43. In Fig. 15.21 the evolution of an initial curve (number 12 
of Fig. 15 . 1 [) under curve shrinking is shown from two points of view. The flow is started at 
r = and stopped at a value of r shortly below T. In Fig. 15.31 the flow of the intersection 
points x int> j(r) corresponding to the initial curves depicted in Fig. 15.11 is shown. 
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Figure 5.2: Plots of the evolution of an N = 1 center-vortex loop (curve 12 of Fig. I5.ip 
under curve shrinking. The thick central line depicts the trajectory of the intersection 
point which coincides with the origin at r = 0. 




Figure 5.3: Flow of the intersection points x intj j(r) for the initial curves depicted in Fig. 15. 1 1 
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5.2.2 Numerical investigation 

In general, the procedure is in close analogy to Sec. 14.2.21 Therefore, we solely mention 
the differences compared to Sec. 14.2.21 The choice of sample points, used to generate the 
initial curves, is now done in such a way that, if interpolated one after another by a cubic 
spline with coinciding initial and end point, the curve will cross itself once. The number of 
grid points n is set to 300 for all curves, except for curve 9 where 500 points are used. For 
the simplicity of the computation of the intersection of the initial curve, the intersection 
point is included in the set of sample points as a double point and is chosen to coincide 
with the origin. 

The search for solutions to the flow equation is proceeded as in Sec. 14.2.21 using the 
Numerical Method of Lines. Starting from r = T\ = the solution is advanced step by 
step as long as r < T or a virtual singularity evolves. If the latter is the case we start 
the same fitting procedure as in Sec. 14.2.21 The error estimate of Sec. 14.2.21 based on the 
constant value of A, is no longer applicable in the case of selfmtersecting curves, however, 
AA(t) = Ai(t) — A2(t) = constant can serve to estimate a value that indicates the validity 
of the numerical solution. The error estimate is defined as 



In the case where both areas of the curve have almost the same value, and the absolute of 
A falls below the numerical precision, we have used 



Since Eq. (|3.27|) involves the two-dimensional curl the straight forward discretization of 
this equation already computes the (signed) difference AA of the areas enclosed by the 
curve. Therefore, we do not have to keep track of the selfmtersection during the numerical 
evolution of the curve. However, the computation of the total area involves knowledge of 
A\ and A2, and this in turn of the intersection point. Once the error estimate exceeds a 
value of 2 the curve is fitted to yield again spatially well separated points. We also used 
the error estimate to recognize the final (non-virtual) singularity at T, where a further 
evolution in the sense of the flow equation is impossible and does not make sense. 

To compute the position of intersection point one searches, at a first step, for those two 
points in the solution set at given r which are spatially nearest to each other, but with 
the restriction that they are element of different line segments of the curve which generate 
the intersection point. Therefore, the spatial distances between the i th and j th point on 
the curve are determined for all pairs of points with \i — j\ > ci m i n , where d m i n depends 
on the considered curve and has to be adjusted for each individually. Then it is searched 
for the least of all distances to find that pair of points which is closest to the intersection 
point. The minimal distance of indices <i m i n is introduced because, as the flow evolves 
the curve, next neighboring points could become (spatially) closer to each other than the 
points nearest to the intersection. Once the pair of points which is next to the intersection 
is found, the two curve segments around these points are approximated by cubic splines. 
Now the intersection of these two splines is computed using Mathematica's FindRoot 

At given r, length and area of the curves are computed with their discrete formulas 
using Mathematica's FiniteDifferenceDerivative. For the same reasons as in Sec. 14.2.21 the 
r-evolution of L, A and x int are interpolated by polynomials. 




(5.8) 



10 6 x (AA( Tl ) - AA(t)). 



(5.9) 
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Figure 5.4: The squares of the coefficients cm(t) entering the ansatz for effective 'action' 



Sm 



A(r) 



1 + 



Cm(t) 

W) 



for T-ordered ensembles up to M = 16. 



Finally, the implicit first-order differential equation dZ/dr = for the coefficient is 
solved using NDSolve for all ensembles sizes and orderings and for both 'actions'. If not 
set at will, the initial condition Co for c(r) was derived from Eq. (|5.7p using Mathematica's 
FindRoot. Variance of mean intersection and entropy were computed. 

A CD-ROM containing the used Mathematica notebooks is attached to the thesi^]. 



5.2.3 Renormalization-group invariance of partition function 

We now present the results of the simulation. For all ensembles Em, the r dependence 
of the coefficient cm in Eq. (15. 3p roughly behaves like a square root oc VTm ~~ r where 
Tm is the weakly ensemble-dependent minimal resolution. For the modified 'action' Sm = 



Lit? 



1 + 



CM ft) 



~Mf) \ x "t A(t) ) ^ ne coefficient cm(t) is well approximated by a linear function oc Tm — t. 
Again, Tm is a weakly ensemble-dependent minimal resolution. For T-ordered ensembles 
the results for cm(t) for the 'actions' Eq. (15.3(1 and Eq. (15. 4h are shown in Fig. 15.41 and 
respectively, in Fig. 15.51 The results for the ensembles E' M do not differ sizably from those 
presented in Fig. 15.41 and respectively, in Fig. 15.51 



5.2.4 Variance of location of selfintersection 

The mean intersection x int (r) over the ensemble Em is defined as 

1 M 

x int (r) = (x int (r),y int (r)) T = — ^ x inM (Y) exp (-S M [^i(r)]) 



(5.10) 



1 Mathematica version 6.0.3 was used. Pay attention to the footnote in Sec. 14.2.21 
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Figure 5.5: The coefficient cm(t) entering the ansatz for the effective 'action' Sm 
' ' 1 + c y/j^ ) for T-ordered ensembles up to M = 16. 



A(r) 



A(r) 



where x intj j(T) — (^int(^")) 2/int(^"))^ is the location of the point of seffintersection of curve 
Xj at r. The scalar statistical deviation Ajv/mt of x int over the ensemble Em is defined as 



AM,kt( T ) = V vaXM,int;a;(T) + var M ,int;j/ (t) 



(5.11) 



where 



var A f, in 



1 M 

— ^2 ( X '^A T ) ~ x int(T)) 2 exp (-S M [xi( T )]) 



i=l 



M 

+ ^E<*W exp(-SM[ Xl (r)]) 



(5.12) 



i=l 



and similarly for the coordinate y. In Fig. 15.61 plots of Ajw jiat (r) are shown when evaluated 
over the ensembles Ei,..., E\§ subject to the 'action' 

s M = m(i + CMiT) 



A{r) 



L(t) 



and the initial condition Luij = 0) = Lm{t = 0). In Fig. 15.71 the according plots of 
AM,mt( r ) are depicted as obtained with the 'action' 



A(t) V" ' Mr) 

and subject to the initial condition Lm{t = 0) = Lm(j = 0). Relaxing the constraint of 
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T-ordering (Em — ► ^Af) does not entail a qualitative change of the results. The fluctuation 
in the first graph of Fig. 15.61 and in Fig. 15.71 representing the trivial ensemble Em=i is 
within the range of the numerical precision. The results presented in Fig. 15.61 and Fig. 15.71 
are unexpected since in the N = sector the variance of the 'center of mass' saturates 
rapidly to finite values. In contrast, for the N = 1 sector, the variance of the location 
of selfintersection initially increases, reaches a maximum, and decreases to zero at a finite 
value of r. This is readily confirmed by the evaluation of the entropy, see next section. 

5.2.5 Evolution of entropy 

Let us now investigate the flow of entropy. The weight-functional Pm is defined as 

Pm{t) = PM[x int ,i(r)] = — exp(-S M [xj(r)]), (5.13) 
and the entropy Em as 

M 

S M (r) = S M [x inM (r)] = ^P A/ [x mtii (r)]log(PM[x mt ,i(r)]) (5.14) 

i=i 

1 M 

= \ og Z M + — VS M [xi(r)] exp(- < S M [x i (r)]) (5.15) 
Zm 

where 5A/[xj(r)] is given by Eq. (|5.1|) . In Figures IB~8l and l5~9l plots are shown for Sm(t) 
when evaluated with the 'action' 



A(t) V L(t) 
and respectively, when evaluated with the modified 'action' 

Sm = l m u + <*m 



for T-ordered ensembles of size M = 1, . . . , 16. The continuous approach of entropy to 
zero at finite values of r implies the spontaneous emergence of order in the system as the 
resolution decreases: starting at a finite value of r, a particular member of Em is singled 
out by its weight approaching unity. This is validated by Fig. 15.101 where the weight- 
functionals Pm are shown for T-ordered ensembles of size M = 2, . . . , 4. The pattern that 
a curve is singled out by its weight-functional as r increases continues for all ensemble 
sizes M . In view of Chapter HJ this behavior is highly unexpected and we conclude that 
the nontrivial topology of the N = 1 sector induces qualitative differences into the coarse- 
graining process. 
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Chapter 6 

Applications 



Here, it should be recalled that a magnetic charge emerging as a result of the dynamically 
broken gauge-symmetry SU(2) — ► U(l) in the deconfining phase is interpreted as an electric 
charge with respect to the U(l)y subgroup of the electroweak sector. In view of Sec. 16.21 
recall that the magnetic center flux of the confining SU(2) e is dually interpreted as electric 
flux. 



6.1 Solitonic fermions 

The notion of a point-like electron has always been plagued by its diverging self-energy. 
Even in classical field theory it is present in the shape of the infinite self-energy of a point 
charge. If the electron is considered to be a little sphere of radius R and mass m e with 
electric charge e attached to the sphere, then the electric field energy U is given by 

U = ^R- <« J > 

Sending R to zero, as we have to if we think of the electron as a point particle, the self- 
energy contribution to the mass of the electron diverges. In Quantum Electrodynamics 
(QED), the problem persists: the correction to the electron mass is still infinite, although 
it has a much softer logarithmic divergence (c = H = 1), 

e 2 

Sm QED = 3^ ln(m e i?). (6.2) 

Therefore, one needs to employ renormalization theory to cope with the emerging diver- 
gences that are to a large extent a direct consequence of locality: the point-particle like 
nature of the electron. and despite the comparative success of renormalisation theory 
the feeling remains that there ought to be a more satisfactory way of doing things." as 
Lewis Ryder put in [8]. 

In the Standard Model, the electron is represented by the famous Dirac equation. 
Though it successfully predicts the electron's antiparticle, the positron, and the magnetic 
moment with an g-factor of 2 it has to introduce the concept of the Dirac sea to make sense 
of the infinite number of negative-energy eigenstates. The Dirac sea leads to an infinite 
contribution to the energy density of the 'vacuum' which has to be canceled, somehow. 
Furthermore, the Standard Model does not provide for a deeper explanation of the value of 
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the magnetic dipole moment other than that following from the Dirac equation and small 
radiative corrections. Moreover, the electron mass enters the QED Lagrangian as a free 
parameter, and the running of which with resolution needs additional experimental input. 

The excitations in the confining phase of SU(2) Yang-Mills thermodynamics are single 
and selfintersecting center-vortex loops. The mass of each intersection point in a self- 
intersecting center-vortex loop is given by the Yang-Mills scale A c . Since a monopole 
(antimonopole) is located at the intersection, it carries one unit of electric charge. Recall 
that a magnetically charged object in the defining gauge theory has to be interpreted as an 
electrically charged object in the Standard Model - and vice versa. In a given segment of a 
flux tube, the monopoles (antimonopoles) can move in both directions: there is a two-fold 
degeneracy of direction of the center flux that is analogue to the two-fold degeneracy of 
the spin projection. Moreover, for each center-vortex loop it is possible to move along the 
entire flux system on a closed curve. Thus, the projection of the dipole moment generated 
by the current of monopoles and antimonopoles inside the vortex core onto a given direc- 
tion in space is two-fold degenerated as well. Therefore, we identify each soliton with a 
spin- 1/2 fermion. Setting the Yang-Mills scale Ac equal to the electron mass m e , we are 
led to interpret N = 1 center- vortex loops as electrons or positrons [TU1 [T6l [23] . 

Let us consider the process of twisting and charge localization more closely. The tran- 
sition from a non-selfintersecting to selfintersecting center-vortex loop is by twisting of a 
non-selfintersecting curve. The emergence of a localized (anti)monopole in the process is 
due to its capture by oppositely directed center fluxes in the core of the intersection (eye 
of the storm). By a rotation of the left half-plane in Fig. I6.1f a) by an angle of 7r, see 
Fig. 16.1T b). each wing of the center- vortex loop forms a closed flux loop by itself, thereby 
introducing equally directed center fluxes at the intersection point. This does not allow 
for an isolation of a single spinning (anti)monopole in the core of the intersection and thus 
is topologically equivalent to the untwisted case Fig. I6,lf a). However, another rotation of 
the left-most half-plane in Fig. 16.11 (c) introduces an intermediate loop which by shrinking 
is capable of isolating a spinning (anti)monopole due to oppositely directed center fluxes. 
Notice that in the last stage of such a shrinking process (short distances between the cores 
of the flux lines), where propagating dual gauge modes are availabl^], there is repulsion 
due to Biot-Savart which needs to be overcome. This necessitates an investment of energy 
manifesting itself in terms of the mass of the isolated (anti)monopole (eye of the storm). 
Alternatively, the emergence of an isolated (anti)monopole is possible by a simple pinching 
of the untwisted curve, again having to overcome local repulsion in the final stage of this 
process. 

In the analysis performed in Chapter El we have solely regarded the situation depicted 
in Fig. 16.11 (d). since the direction of center flux within a given curve segment is irrelevant 
for the process of a spatial coarse-graining microscopically described by the same curve- 
shrinking flow as applied to N = center- vortex loops in Chapter |H 

There are also phenomenological reasons that argue for a non-local nature of the elec- 
tron. Recall that the imaginary part of the pressure in the confining phase starts to 
dominate when approaching the Hagedorn transition, thereby inducing microturbulences 
in the plasma (see Sec. I2.4|) . Such a nonthermal behaviour is likely to be related to the ob- 
served but poorly understood microturbulences and internal transport barriers in tokamak 
experiments with magnetically confined plasma. This presumes to identify the neutrino 



1 On large distances these modes are infinitely massive which is characteristic of the confining phase. 
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9 = 9 = n 9 = n 




(b) (d) 



Figure 6.1: (Topological) transition from the N = sector (a), (b), (c) to the N = 1 sector 
(d) by twisting and subsequent capture of a magnetic (anti)monopole in the core of the 
final intersection. Arrows indicate the direction of center flux. 



and the electron with the non-selfintersecting and the one-fold selfintersecting center- vortex 
loop of the confining phase of SU(2) e with Yang-Mills scale A c = m e = 511 keV. Here, it 
should be noticed that due to the absence of an antiparticle in the case of a N = soliton 
neutrinos need to be of Majorana type which is in compliance with the successful search 
for the neutrinoless double (3 decay [3]. 

This interpretation is also supported by recent high-temperature Z-Pinch experiments 
at Sandia National Laboratories detecting an unexpected powerful contained explosion. 
There, an electric current rises in a wire array up to ~ 20 MA within ~ 100 ns, thereby 
the wire is transformed into a plasma column. The strong magnetic field induced by the 
current results in an inward directed (magnetic) pressure P m which compresses the plasma 
until it collapses. In the course of the implosion, the ions and electrons are accelerated 
towards the plasma axis. The radiated soft x-ray energy is as much as four times the kinetic 
energy that is expected to be released by the intersection of ions and electrons. However, 
before the plasma explodes it stabilizes for about 5 ns (stagnation). The measured electron 
temperature T e is found to be ~ 3 keV at stagnation. Preceding the explosion, an ion 
temperature Tj about 300 keV is sustained shortly after the plasma has stagnated [15J. 
The outward directed plasma pressure P p needs to be equal in magnitude to P m = — 1.8 x 
10~ 12 MeV in order that the implosion stagnates. The measured electron temperature is a 
factor 1/8.5 too low if it is asserted that the plasma pressure is carried by electrons only; this 
would correspond to T e ~ 31.55 keV. In [50], the observed imbalance between the energy in- 
and output was addressed to the rapid (~ 1 ... 2 ns) conversion of magnetic field energy to 
a very-high-ion-temperature plasma by the unexpected forming of short wavelength m = 
magnetohydrodynamics (MHD) instabilities at stagnation which subsequently provide 
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associated viscous ion heating. At stagnation, the ions reach much more rapidly than the 
electrons (~ 1 ... 2 s) a temperature of ~ 300 keV, and subsequently heat the electrons 
up to ~ 300 keV, at least locally By equipartition the ion energy is transferred to the 
electrons, leading to the soft x-ray radiation. 

According to Sec. l2.4l and the discussion in [16] this will involve center- vortex loops with 
a higher number of selfintersections. These accelerate the transit of thermal energy from 
ions to electrons and generate a larger energy density and pressure than expected from 
electron dynamics only. As a consequence, the electron temperature rises rapidly after 
stagnation. After the ion-induced heating the Debye screening mass mp of a conventional 
electron-photon plasma is comparable to T e . Thus, at about ~ 5 ns after stagnation, the 
plasma is absolutely opaque and no radiation is released. 

For T e < 0.6 m e and Nc = 6, the truncated sums pv c = ^2n c =i^c,n and pN c = 
J2n=iPc,n of the pressure P c = J2n=o p c,n and the energy density p c = J2n=o Pc,N 
in the electronic system are in the regime of asymptotic convergence [16J. The case i = 1 
corresponds to a contribution of electrons and positrons only. While for T e <C m e center- 
vortex loops with higher mass (higher number of intersections) are strongly suppressed, 
these do significantly contribute to the pressure and energy density for T e > 0.1 MeV. 
At T e = 0.25 MeV, the relative partial pressure and the relative partial energy density 
is Pq ~ 3 Pi and respectively, pe ~ 4.8 pi, and at T e = 0.3 MeV, already Pq ~ 5 Pi and 
respectively, p& ~ 9.4 p\. Notice that at T e = 0.25 MeV, the ratio of Pi to the magnetic 
pressure P m at stagnation is: jr- ~ —4.4 x 10 8 . The existence of center- vortex loops with 
higher mass facilitates the rapid increase of T e to Tj ~ 0.3 MeV and eventually initiates 
the powerful explosion. 

When the electron temperature approaches a value of about 0.5 MeV the Hagedorn 
transition towards the preconfining phase is expected to take place where all charges con- 
dense densely packed into a new ground state. The Z vector boson of the Standard Model 
is identified with the decoupled dual gauge mode in the magnetic phase of SU(2) e . 

The electron appears to be structureless for (nearly) all external momenta that are 
used to probe the system because of the existence of a Hagedorn-like density of states: 
the invested energy deposited into the vertex is converted into entropy associated with 
the excitations of a large number of unstable and heavy resonances (see Fig. 11.11 for the 
excitations with up to iV = 3 selfintersections). Only for momenta comparable to the 
Yang-Mills scale A c , the BPS monopole located at the intersection becomes excited and 
reveals a part of its structure. For momenta sizeably below A c , there is nothing to be 
excited in BPS monopole. 

6.2 High-temperature superconductivity 

Let us now sketch an alternative approach to high-temperature (high T c ) superconductivity. 
Recall, that the magnetic center flux, dually interpreted as electric center flux, is two-fold 
degenerated. Now the electric charges that travel along the flux lines in the vortex core 
produce a magnetic dipole moment. The projection of which onto a given direction in 
space is either parallel or antiparallel and represents the two-fold degeneracy of the spin 
projection. Here, it should be recalled that a shift of the intersection point of an isolated 
N = 1 center- vortex loop leaves the mass of this soliton invariant. 

Coincidentally, there are quantum systems in nature the unconventional behaviour of 
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which seems to be closely related to the restriction of electron dynamics to two space 
dimensions. In particular, these include high T c superconductors such as the family of 
superconducting materials largely containing (rare-earth) doped cooper-oxide (cuprates) 
planes as well as the recently discovered new class of layered oxypnictide superconductors. 
Let us consider the former at first. In both cases, superconducting layers of magnetic 
moments are interspersed with layers of nonmagnetic material. This nonmagnetic material 
also serves as an reservoir that provides, by doping, for the electrons and screens the 
Coulomb repulsion in the superconducting layer between them. Now the question arises 
how long-range interactions of magnetic moments at given optimal doping and sufficiently 
low temperature lead to superconductivity in the cuprate layers. 

Since, at small enough temperatures, copper-oxide planes are Mott insulators with long- 
range antiferromagnetic order of spins, the conventional Hubbard model must be used. A 
canonical transformation involving a Gutzwiller projection leads to the l t — J' model, where 
t describes the hopping of electrons from site to site and J the superexchange J = 4t 2 /U 
with U describing the Coulomb repulsion. Here, the Gutzwiller projection, which removes 
most of the phonon pairing interaction, is mandatory. Variation of the electronic degrees 
of freedom results in a set of gap equations for the ground state that give the predicted 
d-wave gap and the superconducting order parameter (related to the critical temperature 
T c ) as a function of doping [B"T] . 

Let us now sketch a somewhat speculative approach to high T c superconductivity being 
well aware of our lacking theoretical knowledge on details in this field of research. The key 
idea is already encoded in Fig. 16.11 (d). According to SU(2) Yang-Mills theory, the electron 
represented as a selfintersecting center-vortex loop is a non-local object the magnetic dipole 
moment of which is only loosely related to the localization of its charge: the magnetic mo- 
ment, carried by the vortex core of the flux lines, receives contributions from line segments 
which are spatially far separated (on the scale of the diameter of the intersection) from the 
location of the electric charge. This suggests a system of planar center- vortex loops trapped 
in a two-dimensional layer where the interaction between vortex lines becomes important 
due to an effective screening of the electron charge leading to an ordering effect. In view 
of the reported strong correlations between electrons in two-dimensional superconducting 
systems [52], we imagine a situation as it is depicted in Fig. 16.21 

Due to Ampere's law equally directed electric flux lines attract each other, whereas op- 
positely directed flux lines experience a repulsive force. So, for a given center-vortex loop, 
there is an attractive interaction of four out of six line segments defined by the neighbouring 
electrons while the other two repulse each other. The existence of interactions between flux 
lines that are mediated by the photon is a consequence of the mixing between the gauge 
groups SU(2) e and SU(2) CMB , the latter pertinent to the existence of propagating photons, 
see [TU1 EH]- It should be noticed that the spin projection of a given electron is equally 
directed for two of its neighbours while the other four have oppositely directed spins. This 
supports the observation that high T c superconductivity is an effect not related to s-wave 
pairing [52] ■ An overlap of flux lines would lead to additional intersection points which 
require an extra amount of energy Ac for each intersection and is therefore energetically 
forbidden because the fluctuations in energy density of the system will not allow for the 
creation of an intersection of mass m e = 511 keV. This leads to a repulsive force as the 
spatial distance between adjacent vortex segments vanishes. In order that an equilibrium 
between attraction and repulsion where the intersection point is fixed with respect to its 
neighbours occurs, as it is depicted in Fig. 16.21 one needs a sufficiently low temperature 
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Figure 6.2: The figure, possibly related to the superconducting state in a cuprate, shows 
an array of strongly correlated center-vortex loops tiling the two-dimensional plane. If 
optimal screening of the electron charge located at the intersection point is provided by 
doping such that an attractive interaction between center-vortex loops due to Ampere's 
Law becomes important, then the attractive force between equally directed center flux 
segments could lead to the indicated equilibrium configuration. For a given electron there 
are six neighbouring line segments two of which experience repulsion while the other four 
experience attraction. An overlap of flux lines would create new intersection points, each of 
mass oc m e , which is topologically forbidden, thus leading to a repulsion at short distances. 

(related to resolution) and an optimal screening of the Coulomb repulsion by the surround- 
ing reservoir layers. If the temperature (resolution) falls below a critical value then the 
fluctuations of the intersection points relative to one another will vanish. Applying an 
external electric field parallel to the plane would set the stiff system of locked electrons 
in a collective motion with zero electric resistance. Macroscopically, this situation is illus- 
trated by a stiff table cloth being pulled over the table in a frictionless way. The measured 
pseudo-gap phase in high T c superconductors will be addressed to local distortions in this 
highly ordered state. These distortions require an finite amount of energy and are due to 
insufficient screening and/or to much of a thermal noise. 

Let us now turn to the recently discovered, new class of high T c superconductors that 
are based on oxypnictides (a class of materials including oxygen, an element of the nitrogen 
group (pnictogen), and one or more other elements). These do not seem to exhibit strong 
correlations between the electrons contained in the two-dimensional (FeAs) layers where the 
electron dynamics takes place, see [53]. If the behaviour of a two-dimensional system of non- 
interacting electrons, which are subjected to an environment represented by a parameter 
r, effectively is describable by a coarse-graining process in a statistical ensemble, as it is 
investigated in Chapter 03, then we should likely address the observation that the entropy 
vanishes at finite r to this particular kind of high T c superconductivity. Namely, the 
observation that no variance of the ensemble average of the position of the electric charge 
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is allowed for at a finite resolution is crucial for the statement that the two-dimensional 
system of free quasiparticles is void of any electric resistance. Free quasiparticles in the 
sense that explicit interactions between the electrons in the superconducting layer are 
absent, but the distortions induced by the noise of the environment is fully taken into 
account. Again, r should be a monotonic function of temperature. 



Chapter 7 



Summary 



In this thesis, we have investigated center-vortex loops with and without selflntersection, 
as they emerge in the confining phase of SU(2) Yang-Mills thermodynamics. In a noisy 
environment, center-vortex loops are subject to a spatial coarse-graining due to a motion 
by curvature that is described by a curve shortening flow. In a statistical description of 
ensembles of center- vortex loops which are (locally) embedded into a two-dimensional flat 
plane, we have defined an effective 'action' in purely geometric terms that is governed by a 
renormalization-group flow driven by the curve shrinking. The 'action' possesses a natural 
decomposition into a conformal and a non-conformal factor. 'Observables', such as the 
position of 'center of mass' (JV = 0), or the intersection point (N = 1), are computed as 
ensemble averages of local or non-local operators on the curves. 

We have made the observation that N = center-vortex loops exhibit a second-order 
transition to the conformal limit of vanishing curve length with a critical mean-field ex- 
ponent of the coefficient: on average, center-vortex loops disappear from the spectrum of 
confining SU(2) Yang-Mills theory, thus generating an asymptotic mass gap. The evolution 
of the variance of the initially sharp position of 'center of mass' saturates at finite value 
within a finite decrease of resolution Q, the latter related to the resolving power used to 
probe to the system. These findings bear a strong family resemblance with the unitary 
time evolution of a free particle in quantum field theory. 

Since we believe that N = center-vortex loops play the role of Majorana neutrinos 
[4], the concept of a neutrino rest mass is no longer applicable. Its mass is the result of the 
distortions induced by the environment it is embedded in and depends on the resolution. 
The disappearance of N = center-vortex loops from the excitation spectrum and the 
absence of a corresponding antiparticle would be manifestations of lepton-number violation 
forbidden in the Standard Model of Particle Physics. 

In the case of one-fold intersecting center- vortex loops, we have obtained the unexpected 
result that a statistical ensemble of initial curves evolves into a highly ordered state. That 
is, only a particular member of the ensemble survives the process of two-dimensional spatial 
coarse-graining. As a consequence, the entropy attributed to the ensemble moves to a zero 
value for a sufficient decrease of resolution. 

We have sketched an alternative approach to high-temperature superconductivity based 
on cuprates. The central observation that these depend highly on strong correlations be- 
tween electrons trapped in a flat two-dimensional layer is attributed to an array of bow-tie- 
like simplices (N = 1 center- vortex loops) tiling the plane. We have also speculated that the 
spontaneous emergence of order in an ensemble of planar N = 1 center-vortex loops could 
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be relevant for the recently discovered new class of oxypnictide layered high-temperature 
superconductors that do not seem to exhibit explicit, strong correlations between electrons 
within the superconducting (FeAs) planes. 

In a sense, we have reversed the usage of the renormalization group. Instead of starting 
with a 'physical' action, from which an equation of motion follows, and demanding the 
system to be invariant under renormalization-group transformations, we have defined a 
geometric effective 'action' the coefficient of which is determined by a renormalization- 
group flow driven by a coarse-graining process (curve shrinking). Afterwards, physical 
'observables', such as the center of mass or the localization of electric charge, were computed 
as mean values in a statistical ensemble. In a manner of speaking, we have 'derived' a 
statistically averaged 'equation of motion'. In this context, we regard resolution over time 
(or temperature) as the more fundamental quantity describing a quantum mechanical or 
statistical system. However, this requires the introduction of a model that relates resolution 
to time (or temperature). 

An obvious extension of our statistical approach would be the account of interactions 
between center- vortex loops; especially in view of two-dimensional systems exhibiting high- 
temperature superconductivity. At a first stage, this would include Coulomb interactions 
between the charges of A^ = 1 center- vortex loops localized at the intersection point, and a 
delta-function-like repulsion due the topologically forbidden overlap of center- vortex loops 
(contact interaction). This could be done by adding interaction terms to the effective 
'action' in the partition function that are weighted accordingly. Considering a multitude 
of configurations of initial curves of center- vortex loops, those configurations will be singled 
out the curves of which are most likely to survive the process of coarse-graining. However, 
the nontrivial issue arises how to gain a weight which relates the purely geometric 'action' 
to the one stemming from the electromagnetic interaction between center-vortex loops. 
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